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ABSTRACT 

We report on the cosmological parameters derived from observations with the Cosmic Background 
Imager (CBI), covering 40 square degrees and the multipole range 300 < £ < 3500. The angular scales 
probed by the CBI correspond to structures which cover the mass range from 10 14 M© to 10 17 M Q , 
and the observations reveal, for the first time, the seeds that gave rise to clusters of galaxies. These 
unique, high-resolution observations also show damping in the power spectrum to £ ~ 2000, which we 
interpret as due to the finite width of the photon-baryon decoupling region and the viscosity operating 
at decoupling. Because the observations extend to much higher £ the CBI results provide information 
complementary to that probed by the BOOMERANG, DASI, MAXIMA, and VSA experiments. When 
the CBI observations are used in combination with those from COBE-DMK we find evidence for a 



flat universe, f2t 



1.00 



+0.11 
-0.12 



densities in cold dark matter, fl c dmh 



(1-ct), a power law index of primordial fluctuations, n s — 1.08i ' 10 , and 



0.023lg$$. With the addition 



7+0.10 



0.16±g g|, andbaryons, fl b h 2 
of large scale structure priors the S} c d m h 2 value is sharpened to 0.10l 03 , and we find Q A = 0.67±g 13 . 
In the £ < 1000 overlap region with the BOOMERANG, DASI, MAXIMA, and VSA experiments, the 
agreement between these four experiments is excellent, and we construct optimal power spectra in the 
CBI bands which demonstrate this agreement. We derive cosmological parameters for the combined 
CMB experiments and show that these parameter determinations are stable as we progress from the 
weak priors using only CMB observations and very broad restrictions on cosmic parameters, through 
the addition of information from large scale structure surveys, Hubble parameter determinations 
and Supernova-la results. The combination of these with CMB observations gives a vacuum energy 



estimate of Oa = 0.70 



+0.05 
-0.05' 



a Hubble parameter h = 0.69 ± 0.04 and a cosmological age of 13.7 ± 0.2 



Gyr. As the observations are pushed to higher multipoles no anomalies relative to standard models 
appear, and extremely good consistency is found between the cosmological parameters derived for the 
CBI observations over the range 610 < £ < 2000 and observations at lower I. 
Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

The angular power spectrum of the Cosmic Microwave 
Background (CMB) has emerged as a major arena in 
which our cosmological models and theories of cosmic 
structure formation can be tested. In this paper we esti- 
mate cosmological parameters from observed power spec- 
tra from the Cosmic Background Imager (CBI) and the 
low £ anchor of the COBE-DMK observations, and we 
relate the CBI observations to BOOMERANG, DASI, 
MAXIMA, and earlier CMB experim ents. We have also 
incorporated results from the VSA l)Scott et alJ 120021) 
which were reported at the same time as these CBI re- 
sults. We also use results from large scale structure stud- 
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ies (LSS), supernova observations (SNla), and Hubble 
constant (HST-ft.) measurements to refine estimates of 
key cosmological parameters, both for CBI+DMR alone, 
and for the CBI results in combination with other CMB 
observations. 

During the year 2000 observing season, the CBI 
covered three dee p fields of diameter roughly 1° 
l)Mason et al J 120031 hereafter Paper II), and three mo- 
saic regions, each of size roughly 13 square degrees 
ijPearson et al .1120031 hereafter Paper III). Methods used 
for power spectrum estimation from these in t erfero metry 
observations are described in iMvers et alJ i ffi KKi. here- 
after Paper IV). The results from the 2001 observing 
season, when combined with the observations from 2000, 
will extend the mosaics to 80 square degrees, roughly 
doubling the amount of mosaic data. This will improve 
upon the parameter estimates given here but has yet to 
be analyzed. 

Over the past several years, as the CMB observa- 
tions have improved, the basic Ci features predicted in 
the eighties for inflation-motivated models in which cos- 
mic structure arises from Gaussian-distributed curva- 
ture fluctuations (e.g., Bond & Efstathiou 1987) have 
emerged: (1) A "Sachs- Wolfe" plateau at low multi- 
pole moments, seen by COBE-DMR l|Bennettl 11996ft 
and other large angle experiments. This probes the 



2 



SIEVERS ET AL. 



DAS I 
MAXIMA 




Fig. 1. — Features in the anisotropy spectrum (from Pa- 
per III) . The firs t acoustic peak is seen at high s ensitivity in the 
BOO MERANG INetterfield et alj|2003) . DASI IHalverson et alJ 
1200211 . and MAXIMA ILee et al.l200H) observations, while the sec- 
ond and third acoustic peaks are seen at lower sensitivity (the 
rectangles indicate the 68% confidence intervals on band-power). 
The circles (dark blue) and squares (green) show the odd and even 
binnings of the CBI results from the joint spectrum of the three 
mosaic fields (see Paper III. Note that the two binnings are highly 
correlated with each other and are not independent measurements 
of the power spectrum). The damping tail is clearly seen in the 
CBI spectrum, and, in the region of overlap, all four experiments 
are in excellent agreement, as is discu ssed in § 15.21 The black curve 
is the joint model also discussed in § 15.21 



gravitational potentials at the last scattering surface 
and along the line of sight on scales that were not in 
causal contact at the epoch of photon decoupling at 
redshift z ~ 1100. (2) The long-sought first acous- 
tic peak at £ ~ 200, tentatively first seen by combin- 
ing r esults from a heterogeneo us mix of CMB results 
(e.g.. IBond. Jaffe fc Kno xl 120001) . then in single experi- 
ments. TOCO (|Miller et alJll999l) and BOOMERANG - 
NA North American Test Flight ipauskopf et alJl2000|) . 
This was followed shortly after by the spectacularly de- 
tailed first peak mapping by the BOOME RANG Antarc 
tic fl ight Ijde Bernardis et alJ 120001: iNetterfield et al 
20021) and MAXIMA ijHananv et all 120011: ILee eTaf 



2001); (3) and the detection of the next few peaks 
and dips by BOOMER ANG (Ide Bernardis et all 120021 
Netterfield et all [2002) and DASI ijHalverson et all 
2002TT Observations from the more recent experiments 
are shown in Figure H Following these successes, three 
othe r key ingred ients remained to be demonstrated (see, 
e.g.. lBondlll996l for a review): (4) The peaks and dips 
continue to higher £ at ever diminishing amplitude, in 
a damping tail directly tied to the viscosity in the 
photon-baryon fluid as they decouple, and to the finite 
width of that decoupling region. (5) A necessary and 
highly predictable linear polarization power spectrum, 
fed by the polarization-dependent Thomson scattering 
of the primary CMB anisotropy developed during pho- 
ton break-out from the decoupling region. (6) Secondary 
anisotropies that are an inevitable consequence as waves 
develop nonlinearly, breaking to form collapsed struc- 
tures. 

The CBI power spectra show clear evidence for a de- 
cline consistent with the damping tail (4) , as can be seen 
in Figure ^ The Sunyaev-Zeldovich effect has of course 
been observed in clusters of galaxies at very high sen- 
sitivities, so (6) is there, and there may be evidence of 
a statistical SZE signal from distant clusters in the CBI 



deep field observations (Papers II and VI), but further 
experimental verification is needed. We address some 
issues associated with (4) and (6) in this paper , and 
(6) is treated in more detail in IBond et all ((2003, here- 
after Paper VI). There are many experiments underway 
to address (5) and the first detection of polarization in 
the microwave back ground has re c ently been announced 
iKovac et alJ Ip^HLdTch et all $200% : CBI has been 
observing polarization since September 2002. A longer 
term possibility would provide a 7th pillar, the detection 
of a component attributable to gravity waves in the CMB 
observations. However, unlike the other six items, infla- 
tion models differ substantially in the amount of gravity 
waves predicted, with very small undetectable imprints 
on the CMB being quite feasible. 

One reason that theory can be so definitive is that 
the primary CMB anisotropies probe the linear regime 
of fluctuations and can be calculated in exquisite de- 
tail. The adiabatic inflation-motivated paradigm contin- 
ues to do remarkably well as our knowledge of Cg and 
other parameters improves. Just as powerful though 
is the large reduction in the space of possible theo- 
ries which accurate measurements of the anisotropics 
provide. The positioning of the peaks is a strong 
argument in favor of predominantly curvature (adia- 
batic) fluctuations as opposed to isocurvature ones. 
In addition the multiple peaks are a strong argument 
in favour of coherent (passive) perturbations as op- 
posed to incoherent (active) ones, e.g., those associ- 
ated with co s mic d ef ect theories of structure formation 
iAllen et al.l Il997t lCo ntaldi 1 Hindmarsh. &: Magueiiol 
ll999HTurok. Pen, fc Seliaklll998fl . 

Minimal inflation-based models characterize the pre- 
dictions with a handful of cosmological parameters 
{Sltot, ^A, Qbh 2 , n cdm h 2 , n s , t c , InCio}. These 
are now very familiar to a strophysicists (see , e.g., 
Bond. Efstathiou & Tegmarkl 119971: Uaffe et all 120011: 



Lange et al.l 1200 It INetterfield et all 120021: iPrvke et all 



2002). The present-day density pj of a component j is 
Slj = 8TrGpj/(3H$) where H„ = lOO/i km s" 1 Mpc" 1 is 
the Hubble constant. We use the notation ujj = Ujh 2 for 
the related physical density, which is more relevant for 
the CMB. We therefore consider the densities of baryons, 
uJb, and of cold dark matter (CDM), w c d m , with total 
matter density parameter £l m = Q c dm + ^fc- The energy 
density parameter associated with a cosmological con- 
stant is fl\ and the total density is f2 to t = + ^A, 
related to the curvature energy density parameter by 
r2fc = 1 — f^tot- The initial spectrum of density perturba- 
tions is described by an amplitude, which we usually take 
to be a factor C\q multiplying the CMB spectrum, and 
the spectral tilt of scalar (density) perturbations, n s (de- 
fined so that the initial three-dimensional perturbation 
power in gravitational potential fluctuations per In k is 
oc k n " , with n s = 1 thus giving a scale invariant spec- 
trum). We also consider the Thomson depth tc to the 
epoch of reionization, presumably associated with lumi- 
nous star formation in the earliest forming dwarf galax- 
ies. Many more parameters may be needed to completely 
describe inflationary models. These include the gravity- 
wave induced tensor amplitude and tilt, variations of 
tilt with wavenumber, relativistic particle densities, more 
complex dynamics associated with the dark energy 
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etc. We discuss these briefly below, but concentrate on 
the minimal set of 7 parameters for this paper. The grid 
of parameter values we have used i n this work is an ex - 
tended version of the one used bv lLange et alJ ([2001), 
and is reproduced in Table 

The main target of the CBI, as with most other exper- 
iments to date, is to measure primary anisotropies of the 
CMB, those which can be calculated using linear pertur- 
bation theory. The maps shown in Papers III and VI are, 
to a first approximation, images of damped sound wave 
patterns that existed about 400,000 years after the Big 
Bang at a time when the photons were freed from the 
plasma. However, the images are actually a projected 
mixture of dominant and subdominant physical processes 
occurring through the photon decoupling "surface" , a 
fuzzy wall at redshift z^ cc ~ 1050, when the Universe 
passed from optically thick to thin to Thomson scattering 

— 1/2 

over a comoving distance ~ 7 ui m Mpc. (Specific num- 
bers here are appropriate for a ACDM universe preferred 
by the parameter estimates in this paper; uj m ~ 0.14 is 
one of our conclusions.) Prior to this epoch, acoustic 
wave patterns in the tightly-coupled photon-baryon fluid 
on scales below the comoving "sound crossing distance" 
at decoupling, ~ 50 w™ 1 ^ 2 Mpc (i.e., ~ 50w,7i 1 ^ 2 kpc 
physical), were viscously damped and strongly so on 

— 1/2 

scales below the damping scale ~ 4cj m Mpc. This is 
closely related to the thickness over which decoupling oc- 
curred. Subsequently, the photons freely-streamed along 
geodesies to us, mapping (through the angular diameter 
distance relation) the post-decoupling spatial structures 
in the temperature to the angular patterns we observe 
now as the primary CMB anisotropies. For example, the 
sound crossing and damping scales translate to multi- 
poles ~ 110 and ~ 1300, respectively. Free-streaming 
along our (linearly perturbed) past light cone leaves the 
pattern largely unaffected, except for the effect of tem- 
poral evolution in the gravitational potential wells as the 
photons propagate through them which leaves a further 
AT imprint, known as the integrated Sachs- Wolfe effect. 

A number of effects complicate this simple picture of 
direct mapping of acoustic compression and rarefaction 
regions: anisotropies are also fed by the electron flow 
at decoupling leading to a Doppler effect, the chang- 
ing gravitational potential as the universe passes from 
domination by relativistic to nonrelativistic species and 
small terms associated with polarization. The damping 
is of course a major radiative transfer problem, connect- 
ing the tightly-coupled baryon-photon fluid regime when 
shear viscosity and thermal conduction accurately de- 
scribe the damping through to the free-streaming regime 
when full transport is needed. Intense theoretical work 
over three decades has put accurate calculations of this 
linear cosmological radiative transfer problem on a firm 
footing. We discuss this further in § 

Of course there are a number of nonlinear effects 
that are also present in the maps. These secondary 
anisotropies include weak-lensing by intervening mass, 
Thomson-scattering by the nonlinear flowing gas once 
it became "reionized" at z ~ 10-20, the thermal and 
kinematic SZ effects, and the red-shifted emission from 
dusty galaxies. They all leave non-Gaussian imprints 
on the CMB sky. Theoretical predictions based on the 
best-fit models suggest that for CBI, only the thermal 



SZ effect would be within striking distance. This effect 
is addressed in Paper VI. 

The structure of the paper is as follows. In § we 
first summarize the various priors from non-CMB obser- 
vations that we use when estimating cosmological param- 
eters. In § Eland Appendix lAl we describe the methods 
we have used for parameter determination from power 
spectrum estimates, and the tests we have performed 
to demonstrate accuracy. We obtain estimates for our 
minimal inflation-based parameter set in § and also 
apply the techniques of Appendix ^ to determine opti- 
mal bandpower spectra in § El We check for consistency 
between the mosaic and deep field results by comparing 
optimally calculated spectra for various subsets of the 
data. In the parameter estimations we restrict ourselves 
to the 1 < 2000 region. The spectrum at £ > 2000 is dis- 
cussed in Paper VI together with the level of secondary 
anisotropy expected from the Sunyaev-Zeldovich effect in 
that regime, and whether this is the origin of the excess 
power seen in the CBI observations at high I. In § El 
we compare our results with those of other CMB exper- 
iments as a further consistency check. We then combine 
all available results with the CBI observations to obtain 
an all inclusive set of parameter estimates in § El In § 
we summarize the physical effects probed directly by the 
CBI observations. Our conclusions are presented in § El 

2. PRIOR PROBABILITIES USED IN COSMOLOGICAL 
PARAMETER EXTRACTION 

We apply a sequence of increasingly strong "prior" 
probabilities successively to the likelihood functions. 
These are well known from the CMB literature (e.g., 
Bond fc ,Iaff3 Il998t I.Taffe et al.ll200l lLange et al J 120011 
Netterfield et al J 12002(1 . The priors we apply are listed 
below. 

1) The "weak-/i" prior: this prior restricts the Hubble 
parameter to 0.45 < h < 0.9, and also imposes an age 
restriction, to > 10 Gyr, and a restriction on the matter 
density, fi m > 0.1. These are all weak priors that most 
cosmologists would readily agree on. 

2) The "flat" prior: as we shall see, the CBI observa- 
tions strongly support Qk ~ 0, so the addition of a flat 
prior to weak-h seems reasonable. This is especially so if 
the target is parameters associated with inflation mod- 
els. Although it is possible for inflation models to give 
large mean curvature with non-negligible |f2fc|, they are 
rather baroque. 

3) The "LSS" prior: The LSS prior we use here is s l ightly 
modified over that used earli er l|Bond fc J affel 119981 
.laffe et all EfflOlt lLange et all I200H iNetterfield et all 
2002) and is described in detail in Paper VI. It involves a 
constraint on the amplitude erf and shape T c g of the (lin- 
ear) density power spectrum. Here as is the rms density 
power on scales corresponding to rich clusters of galaxies 
(8/i _1 Mpc) and r e ff mainly parameterizes the critical 
length scale when the universe passed from dominance 
by relativistic matter to dominance by non-relativistic 
matter. Both constraints depend upon our basic mini- 
mal parameter set in complex ways. The distributions in 
both these LSS parameters are taken to be quite broad, 
akin to a "weak LSS prior". We take a distribution for 
the combination ogfi^ which is a Gaussian (first error) 
smeared by a uniform (top-hat) distribution (second er- 
ror): 0.47±;Sl;jZ;o|. The 0.47 value is about 15% below 
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Table 1. Parameter Grid for Likelihood Analysis 



Parameter Grid: 





0.9 


0.7 


0.5 


0.3 


0.2 


0.15 


0.1 


0.05 





-0.05 


-0.1 


-0.15 




-0.2 


-0.3 


-0.5 




















u cdm 


0.03 


0.06 


0.08 


0.10 


0.12 


0.14 


0.17 


0.22 


0.27 


0.33 


0.40 


0.55 




0.003125 


0.00625 


0.0125 


0.0175 


0.020 


0.0225 


0.025 


0.030 


0.035 


0.04 


0.05 


0.075 




0.10 


0.15 


0.2 




















n A 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


n s 


1.5 


1.45 


1.4 


1.35 


1.3 


1.25 


1.2 


1.175 


1.15 


1.125 


1.1 






1.075 


1.05 


1.025 


1.0 


0.975 


0.95 


0.925 


0.9 


0.875 


0.85 


0.825 






0.8 


0.775 


0.75 


0.725 


0.7 


0.65 


0.6 


0.55 


0.5 








Tc 





0.025 


0.05 


0.075 


0.1 


0.15 


0.2 


0.3 


0.4 


0.5 


0.7 





the value adopted in the earlier CMB papers, as discussed 
in Paper VI, but the distribution shape is the same. To 
use it, the relation between Cio and a\ is needed for each 
cosmological model considered. For the shape prior, we 
make use of the similarity, over the wavenumber band 
that most large scale structure data probes, between 
changes in the spectral index n s and changes in the 
shape parameter T w fl m h exp[— Os(l + il ln 1 (2/i) 1 / 2 )]. 
The latter includes a strong dependence upon £l m h as 
well as a rough Ob modification. We combine the two 
dependences into a single constraint on the parame- 
ter r off = r + (n s - i)/2 with r eff = o.2it;8S;t:2S. 

a broad distribution over the 0.1 to 0.3 range. It is 
slightly less skewed to lower r c g than the distribution 
used in the earlier studies. These small changes, which 
acco rd better with the emergi ng LSS data from the 
2dF peacock et all2001l). SPSS ISzalav et al.l200lD and 
weak lensing survey s ( Baconetjdjj2j)02jJB^ 
120021 IRefregier etaLl 12002 Ivan Waerbeke et all 1200^ 7 
have a very small impact on the cosmic parameters we 
derive when the LSS prior is applied, mainly because the 
CMB results are now so good that the LSS prior is not 
as powerful a delimiter as it used to be. 

4) The "HST-h" prior: the HST key project has led to 
more restrictive estimates of the Hubble parameter us- 
ing Cepheid data: h = 0.72±0.08 ijFreedman et all2001t 
iMould et al. 2000), where these are Gaussian 1-sigma er- 
rors. We denote this in the paper by the HST-/i prior. 

5) The "SN" prior: Comparison of observations of a large 
number of distant and nearby supernovae of Type la 
leads to a constraint i n the £l m —Q a plane, or equivalently 
the Q k — f^A plane ijPerlmutter et all Il999t IRiess et al.l 
1998), independent of the other cosmic parameters in 
our minimal set. 

3. PARAMETERIZED POWER SPECTRA FROM RADICALLY 
COMPRESSED BANDPOWERS 

We wish to determine likelihood functions C(y a ) = 
P(D\y a ) for data sets D as a function of parameters y a . 
These can be cosmological, as in the minimal inflation 
set described above; bandpowers for a discrete binning 
of Ci spectra; or experimental, as for calibration and 
beam uncertainties. We consider two classes of param- 
eters, those which are constrained at prescribed values 
(external) and those which we allow to dynamically re- 
lax to their maximum likelihood values (internal). For 
the cosmological parameter set, we treat all but Cio as 
external, with h\C(y a ) determined on a 8-million point 



grid. Calibrations, beam uncertainties, bandpowers and 
Cio are treated as internal, with error estimates in the 
neighbourhood of the maximum likelihood made from 
the curvature (second derivative) matrix evaluated there. 
For example, for the optimal bandpower case we treat in 
§ |5J the parameterization is of the simple form: 

c t = Y,y h c[ s) ^M (i) 

b 

(s) 

in terms of shapes C) and window functions ipbi which 
define a partition of unity (X^V'M = !)■ An obvious 
choice for tpbe is the top hat \u, defined to be unity for 
lb < ( < (b+i and zero outside. It is indeed the one we 
use for optimal spectra. 

Ideally we would use all of the information available, 
e.g., a pixel map with errors described by a pixel- pixel 
correlation matrix, to compute C This has been done in 
the past for limited parameter sets, but the parameter 
spaces we treat now are large enough that a large algo- 
rithmic speedup would be needed. In practice, we first 
go through a stage of "radically compressing" the infor- 
mation into a set of bandpowers with errors. However it 
is essential for accuracy that the entire likelihood surface 
be well represented. We show that this is true for the 
CBI observations. Our data analysis pipeline, described 
in Paper IV, grids the visibility data into estimators, the 
covariances of which are also computed, and then de- 
termines the maximum likelihood power spectrum and 
the curvature of the likelihood function about that max- 
imum. In the last step, the power spectrum for a given 
experiment is parameterized by a discrete sum over con- 
tiguous bands B as in equation (JTJ, with q B replacing 
y b and \Bt for tj)u- For the power spectra estimates of 
Papers II and III, a flat shape (C^f 1 constant) was used. 

Even though {xBe} may be chosen for the input win- 
dow functions in the parameterized model of equation 
processing through the actual (u, v) coverage results 
in an effective set of window functions {Ws(l)} (Pa- 
per IV). These are plotted for the deep and mosaic ob- 
servations in Papers II and III. The Wb {&) depend upon 
signal-to-noise ratio as a function of i and spill over to i 
values that lie beyond the support of the top hat xbi- 

The result of the "radical compression" of the full CBI 
noisy visibility data set, as described in Paper IV, is: 
{q B ,qSo ifes, 9n. {F'^BB', VBi}- Here {q B } are max- 
imum likelihood values of {q B }, (F~ 1 )bb i is the in- 
verse Fisher matrix which would describe the correla- 
tions among the bandpowers if the likelihood distribu- 
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tion were Gaussian, and ipei are the bandpower win- 
dow functions that convert a spectrum Ct to expected 
bandpower values {q B }- In addition to the bandpowers, 
there are often extra parameters associated with other 
contributions to the signal that must be simultaneously 
determined. For example, with the CBI data, we have 
two classes of sources which we take into account: NVSS 
sources having flux densities greater than 3.4 mJy at 1.4 
GHz, and "residual" faint sources which we model as an 
isotropic Gaussian random field. The NVSS sources have 
known positions, which enable us to define point-source 
template structures on the data. Rather than trying to 
model the source amplitudes in detail, we simply project 
them out of the data by making the multipliers, g src , 
of the templates very large. We estimate the amplitude 
of the residual source contribution by extrapolating the 
CBI 31 GHz source counts to fainter flux densities. 

We usually fix the uncertainty in those amplitudes as 
well, but in some tests have allowed it to vary. In that 
case, a nuisance parameter is introduced which increases 
the Fisher matrix dimension by one. Forming (F~ 1 )bb > 
and using it in the treatment of the data is equivalent 
to having marginalized over the nuisance parameter, so 
in effect (F~ 1 )bb i is what is needed in all cases. If we 
target a limited number of bandpowers for parameter es- 
timation, {F~ 1 )bb i is truncated to those bands. This is 
mathematically identical to treating all bandpowers, but 
marginalizing (integrating) over the "unobserved" bands 
we have cut out. 

The effective or generalized total noise in each band, 
?Nt = 9n + <?ras> includes the noise itself q B and the 
source contributions, q B c and q B s . Estimation of these 
noise and source bandpowers is performed after the max- 
imum likelihood bandpowers are found, and they are cal- 
culated within the same calculational framework, as de- 
scribed in Paper IV. To calculate theoretical bandpowers 
from a given Ci, we also need to specify a set of window 
functions relating the contribution of multipole £ to band 

B, ifBl- 

The offset lognormal approximation we use to charac- 
terize the likelihood surfaces is described in Appendix lAl 
and shown to be an excellent approximation to our CBI 
likelihood functions. 

The CBI £ range and binnings we have used are given 
in § 0] The other experimental results we use have a 
variety of £ coverage and binning: BOOMERANG cov- 
ers an I range of 75-1125 with band width A£ = 50, 
DASI covers 104-864 with A£ variable between 70 and 
100, MAXIMA covers 73-1161 with width 75, and DMR 
covers low £, from 2 to 30 (although we start from 3 
because the quadrupole has a Galactic contamination). 
Bands for TOCO, BOOMERANG-NA, and VSA are de - 
scri bed bv iMiller et all lfl999h . iMauskopf et"aD |2000), 
and lScott et all l|2002f) . The "Apr99" combinatio n of ex - 
periments was introduced in lBond. Jaffe fc Knoxl |2000). 
Together with CBI mosaic, these make up the "all-data" 
combination which we use extensively when combining 
CBI with other experiments. 

The results that we have used, and the labels by which 
we refer to them, are summarized in Tabic [3 

4. COSMOLOGY WITH THE CBI AND ROBUSTNESS TESTS 

In this section we present the cosmological results we 
have derived from the CBI observations, and compare 



these results with those from non-CMB observations. In 
§ El wc compare the CBI observations and parameters 
with those derived from other CMB experiments and we 
combine all of these CMB data to determine the best 
overall values of the cosmological parameters that can be 
derived from the CMB data when combined with large 
scale structure studies, the HST H a project, and super- 
nova type la results. 

In the derivation of cosmological parameters from the 
CBI observations we have carried out a large number 
of consistency tests and checks. We find excellent con- 
sistency in all of the tests we have performed. Some of 
these tests are described in S I4.1I and the more important 
of these tests are summarized in § 14.21 

4.1. The Primary CBI Results 

The basic set of 7 parameters for our fiducial mini- 
mal inflation model, {fitot, &bh 2 , O c d, n /i 2 , n s , Tc, 
InCio}, is described in the introduction, and the grid of 
these parameters is given in Tabled The amplitude C 10 
is a continuous variable. The effect on parameter de- 
terminations of the database boundary and the various 
priors de fined on the space that we use is described in 
detail bv lLange et"aTl l|200lD . 

As described in Paper III, our standard parameter 
determinations for the mosaic observations have been 
made with bins of width A£ = 140, with two alter- 
nate locations of the bins. The "even" binning has 
£ B = 260+ 140B (1 < B < 23), while the "odd" binning 
has £ B = 190 + 1405 (1 < B < 23), where £ B is the 
upper limit of the bin. Data derived with these binnings 
are denoted by CBIel40 and CBIol40 in this paper. A 
coarser binning used has width A£ = 200, with "even" 
spacing £ B = 200 + 2005 (1 < B < 16) and "odd" spac- 
ing £ B = 100 + 200B (1 < B < 16), with the correspond- 
ing data denoted by CBIe200 and CBIo200. Correla- 
tions are strongest between adjacent bins and are typ- 
ically negative for interferometry data. The maximum 
anti-correlation between adjacent bins is about 25% for 
A£ = 140 and about 15% for A£ = 200. The CBIdeep 
standard bins begin at 500, 880, 1445, 2010, 2388, 3000, 
ending at 4000. There is a lower £ bin as well, but we 
do not include it in parameter analysis (we marginalize 
over it). 

The primary results of the CBI cosmological param- 
eter extraction are shown in Table 01 and Figure El 
These show parameter determinations after marginaliza- 
tion over all other parameters for the CBIol40+DMR 
dataset and for all combinations of the priors described 
in § EI 

4.1.1. The Geometry and the Primordial Fluctuation 

Spectrum 

We begin the cosmological parameters discussion by 
considering 17 t ot and n s . We see from Table El that 
under the weak-h prior assumption the combination 
CBIol40+DMR yields Q tot = 1.00 and n s = 

To illustrate clearly the constraints that the CBI ob- 
servations impose we show in FigureElthe likelihoods ob- 
tained from DMR alone under the various priors. Note 
that the fik likelihood curve is very broad. This figure 
shows that the above tight cosmological constraints do 
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Table 2. Data Sets 



Label 



Data Set 



All-data 
Apr99 

BOOMERANG-NA 

BOOMERANG 

CBIol40 

CBIc-200 

CBIel40 

CBIe200 

CBIol40(£ > 610) 

CBIdeep 

DASI 

MAXIMA 

TOCO 

VSA 



Apr99, BOOMERANG-NA, BOOMERANG, CBIol40, DASI, DMR, MAXIMA, TOCO, VSA 
Compilation of 17 experiments prior to April 1999, by JBonci^Jaffe^^Knox (2000J 
North American Test Flight of BOOMERANGfMau skopf et al. 2000*1 

Antarctic Flight of BOOMER ANG Cde Bernardis et al. 2000. 2002: Netterfield et al. 2002: R.uhl et al. 
CBI 02h+14h+20h mosaics, odd bins, A£ = 140 
■20h mosaics, odd bins, A£ = 200 
20h mosaics, even bins, A£ = 140 
-20h mosaics, even bins, A£ = 200 

-20h mosaics, odd bins, A£ = 140, band-powers at £ < 610 discarded 



2002) 



CBI 02h 
CBI 02h 
CBI 02h 
CBI 02h 
CBI 08h 



-14h- 
-14h- 
-14h- 
-14h- 
-14h- 



-20h deep fields 



fHalverson et al. 2002'! . 
^^^^^^^^^^001; Lee et al. 
(Miller et aTmfGT ' 
(Scott et al"*2002) 



2001) 



Note. 



Labels used in the text to designate different data sets used in the data analysis and data tests 



Table 3. Cosmic Parameters for Various Priors using CBIol40+DMR 



Priors 
wk-/i 

wk-h+LSS 

wk-/i+SN 

wk-Zi+LSS+SN 



1.00 

i.05: 

1.03| 
1.04! 



1 07"l" 
' 8:11 



HQ 1 ') 
0.029»;M| 

0.029» : «« 



x h 2 n A 




0.58! 
0.60! 
0.65! 
0.64! 



8:11 
Mi 
8:1? 



Age 

13.9H 
15.4| : | 
14 7^4 
15.0?'? 



< 0.66 

< 0.66 

< 0.67 

< 0.67 



rlat+wk-/i (1.00) 

rlat+wk-/i+LSS (1.00) 

flat+wk-/i+SN (1.00) 

flat+wk-h+LSS+SN (1.00) 



1 07 - 11 
1.05°^ 

1.10: 

1.08! 



8:11 



0.023°;°!° 0.15° 
0-024«;MI ou n 
0.025°'Sl O.llS 



0.025 



n.oof) 



0.11 



A7 - 25 
Q-27 
0.67° lc 

7C&'-$ 
Q-Q7 
710 06 
u -' 1 0.06 



0.54 



,,0.24 
1.24 
34 u l 2 
QJ.2 

0.30RR; 



0.29! 



O.Ofi 



0.068°;|| 

u - uo '0.O2Q 
n nzoO-Olt 
Q017 

n ncou.016 



0.60 



i0.12 
(I 66^ : " 
7()" M 



0.69 



().()<) 



14.0} ! < 0.65 

14.2]| < 0.62 

13.8| : | < 0.65 

13.9}| < 0.63 



flat+HST-/i (1.00) 

flat+HST-h+LSS (1.00) 

flat+HST-h+SN (1.00) 

flat+HST-h+LSS+SN (1.00) 



1.09°;" 0.0262 

im°:ll 0.026° 

0.0268 
0.0268 



i.io; 
Los; 



8:11 



0.13° 
0.11» 
0.12° 

0-U8 



0.658 
0.71° 
0.71° 
0.7lg 



0.38! 
0.29! 
0.29| 
0.29; 



0.0582 
0.0542 
0.0522 
0.0542 



81? 

017 



0.1 

0.70; 
0.71; 

0.70! 



13.31? < 0.65 
,1:? 



13.8J, 
13.6H 
7 l:i 



< 0.64 

< 0.65 
13.7} } < 0.63 



Note. — Estimates of the 6 external cosmological parameters that characterize our fiducial minimal-inflation model set as progressively more 
restrictive prior probabilities arc imposed, (re is put at the end because it is relatively poorly constrained, even with the priors.) Central values and 
lcr limits for the 6 parameters are found from the 16%, 50% and 84% integrals of the marginalized likelihood. For the other "derived" parameters 
listed, the values are means and variances of the variables calculated over the full probability distribution, wk-h requires 0.45 < h < 0.90, Age 
> 10 Gyr, and Q m > 0.1. The sequence shows what happens when LSS, SN and LSS+SN priors are imposed. While the first four rows allow Qtot 
to be free, the next four have Otot pegged to unity, a number strongly suggested by the CMB data. The final 4 rows show the "strong-V prior, a 
Gaussian centered on h — 0.71 with dispersion ±0.076, obtained for the Hubble key project. When the lcr errors are large it is usual that there is 
a poor detection, and sometimes there can be multiple peaks in the ID projected likelihood. 



not arise from the priors. Rather, the other experiments 
nicely complement the CMB, greatly enhancing the dis- 
criminatory power of any single data set. These likeli- 
hood curves should also be contras ted with the "prior - 
only" likelihood curves presented bv lLange et alJ |2001). 
The results of adding the CBI data are shown in Figure[3 
The effect of combining the CBI data with the DMR data 
under the LSS prior is to reduce the uncertainties in f2 tot 
— with these priors f2 to t = l-05^QQg. Note that the un- 
certainties on n s are almost independent of the priors, 
and are always ~ 10%. 

The CBI has very little sensitivity at £ < 400, thus 
these determinations of fitot and n s are basically inde- 
pendent of the first acoustic peak. We have explored the 
degree to which the CBI results depend on the low-£ data 
by eliminating the band-powers at £ < 610, and running 
the same sequence of tests of increasingly restrictive pri- 



ors. The results for the cut at £ = 610 are shown in 
Figure 0] and in Table The likelihood is still sharply 
peaked at f2 to t ~ 1 even though the data from the region 
of the first and second peaks has been discarded. 

The effect of discarding the data at £ < 610 on n s 
is seen in Figure 01 We see that the likelihood again 
peaks near unity, showing that the determination of a 
near scale-invariant fluctuation spectrum is not depen- 
dent on the first or second acoustic peaks in the CBI 
data. 

We conclude that both the geometry and the fluctua- 
tion spectrum are highly consistent with the predictions 
of the minimal inflationary theory, and that this con- 
sistency applies even when the data at £ < 610, corre- 
sponding to the regions of the first and second acoustic 
peaks, are discarded. These results on fi t ot and n s are 
therefore independent of previous results based on ob- 
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10 0.2 0.4 




Fig. 2. — ID projected likelihood functions calculated for the 
CBIol40+DMR data. All panels include the weak-h (solid dark 
blue) and LSS+weak-h(short-dash-dotted red) priors. (LSS is 
the large-scale structure prior.) The Qj, panel also shows what 
the whole Q-database gives before the weak-h prior is imposed 
(black dotted). We note that even in the absence of CMB data 
there is a bias towards the closed models lLangc ct al. 2001). In 
the other panels, flat+weak-h (long-dashed-dotted light blue) and 
LSS+flat+weak-h (dashed green) are plotted. Notice how sta- 
ble the n s determination is, independent of priors. We see here 
that, under priors ranging from the weak-h prior to the weak- 
/i+LSS+fiat priors, the CBI provides a useful measure of four out 
of the six fundamental parameters shown. This is independent of 
the first acoustic peak, where the CBI has low sensitivity, and is 
also largely independent of the spectrum below I ~ 610 for all but 
f!;,/i 2 (see text). 



1 

, 0.8 
N 8 0.6 
^ 0.4 
0.2 




:■. /A whole : 

r'"if}\ wk -j 

L // LSS+wt 
- if i \ \ 

7 \ \ ' ~ 


' ' • \ ' ' • • \- 

-p/ - 

r_.fiaHfjfc A -. 

- LS$4flat+4\ ~ 




~ \ ^ 



-0.5 0.5 



0.5 



10 0.2 0.4 
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n h h 2 n„h 2 n . 



Fig. 4. — Same as Fig. El except that CBI data at I < 610 have 
been discarded. We see that, under the weak-h and LSS+weak-?t 
priors, n k is peaked near zero, as is the case for the whole data 
set, showing that in the CBI data the evidence for a flat universe 
is coming from the Grange above the second acoustic peak where 
damping dominates, as well as from the lower-^ range. We see here 
that under increasingly restrictive priors the CBI data at I > 610 
provide useful constraints on n a , £l m h 2 and Qai showing that these 
results are driven by the shape (and level) of the spectrum at high- 
l independent of the results at low-^. Note, however, that when 
the CBl data are restricted to i > 610 they do not provide a useful 



measure of H^/i 2 (compare with Fig. |2*1, 




0.02 0.04 0.1 0.2 0.3 0.5 1 

n k h 2 n.h ! n . 



Fig. 3. — Results obtained using DMR alone. This gives an idea 
of the role of the LSS prior in sharpening up detections for DMR. 
Note that DMR did reasonably well by itself in firs t indicating for 
this class of models that n s ~ 1 fe.g.. lB*ondlll996T) . Of course it 
could not determine u)^ and the structure in Qfc an d C2a can be 
traced to Cg -database constraints (Langc ct al. 20(Mj). Comparison 
with Fig. |21 shows the greatly improved constraints when the CBI 
data are added. 



servations over the multipole range of the first and sec- 
ond acoustic peaks. The I < 610 and I > 610 regions 
of the angular spectrum indicate a consistent power law 
for the primordial spectrum, for the minimal models we 



consider, over the range of scales covered by the obser- 
vations which now extend down to the scales of clusters 
of galaxies probed directly by LSS studies. 

4.1.2. The Non-Baryonic and Baryonic Matter 
Densities and the Cosmological Constant 

The constraints on the density in dark matter, w c d m , 
and the cosmological constant, f^A are shown in Fig- 
ure [3 One sees that these are tightly constrained when 
the LSS prior is added. The effect of adding the CBI 
data can be seen by comparing the weak-/i+LSS prior 
results in Figure El with those in Figure El We see that 
the CBI data reduce the uncertainties dramatically, and 
that u cdm = 0. 101^3 and n A = 0.67l^;^ for the weak- 
ft.+LSS priors case. If, in addition, we assume a flat 
geometry, the uncertainties in the dark matter density 
are further improved, w c d m = O.IIIq'qji while those for 
the cosmological constant are the same. We see from 
Figure 01 and from Table 01 that the constraints on the 
dark matter density are also reasonably tight when the 
data at i < 610 are discarded. The results on J7a are lit- 
tle changed when the data below £ = 610 are discarded 
although the uncertainties are larger, as can be seen by 
comparing Figure 01 with Figure |21 Thus these results on 
Wcdm and Oa are also established over the high-£ range 
independent of the first and second acoustic peaks. 

The fractional CBI constraints on the baryonic matter 
density are not nearly as tight as those on non-baryonic 
matter. With the weak-Zi+LSS+fiat priors, we see from 
Table 01 that u b = 0.024±™H. This is co nsistent with 
the results from Big Bang nucleosynthes is ijBurles et alJ 
l2000llOlive et al.ll99fllTvtler et al.l2 000). Whenwedis- 
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card the data below I = 610, we do not yet have a use- 
ful constraint on f2f,/i 2 (see Figure UJ. Although weaker, 
the potential constraints on the baryonic matter density 
from the high-€ region of the spectrum are important be- 
cause this is, in principle, a direct and independent way 
of measuring flbh 2 . 

From the combination of the densities in both baryonic 
and non-baryonic matter, we see that the CBI provides 
compelling evidence for a matter density significantly 
lower than the critical density required to close the uni- 
verse, with n m = 0.34 ±0.12, with n b = 0.057 ±0.020 in 
baryons. This result, combined with the flat geometry, 
requires a significant energy component of the universe 
to be supplied by something other than matter, and this 
we assume under our minimal inflationary scenario to be 
the cosmological constant. 

4.1.3. The Hubble Constant and the Age of the 

Universe 

Significant measures of the Hubble constant and the 
age of the universe are again obtained under the weak- 
/i+LSS+flat priors (see Table 0). We find, under these 
priors, that h = 0.66 ± 0.11 and t = 14.2 ± 1.3 Gyr. 
These are in excellent agreemen t with recent determi- 
nations of the Hu bble constant (Frccdr nan et alJ 120011: 
iMould et alJ l2000]) and the ages of the oldest stars in 
globular clusters. 

4.2. Further Robustness Tests 

We have carried out a large number of tests, further to 
those described in the previous section, of the parameter 
extraction from the CBI observations. We describe some 
of the more important of these tests in this section. 

The effects of discarding the CBI data at I < 610 on 
the full suite of parameters and priors can be seen by 
comparing Tableland Tabled In this cut, we discarded 
the first three bins of CBI data. It can be seen here 
that the constraints on cosmological parameters degrade 
gracefully as data are discarded. We have also tested 
the effect of discarding the first four bins of CBI data 
(i.e., up to £ — 750) and we find that the uncertainties 
continue to increase, as expected, but no large variations 
in the central values of the parameters are seen. It is 
clear, therefore, that the CBI results are robust in this 
regard. 

We have compared the results derived from the two 
alternate binnings of the data — the "odd" binning and 
the "even" binning. The results are shown in Figure [5] 
Figured an d Tabled We see that the derived values 
of the parameters agree to within the uncertainties in all 
cases, and that the uncertainties are comparable. This 
comparison demonstrates clearly that there is no depen- 
dence of the cosmological results on the binning choice. 

We have also compared the values of the parameters 
extracted using bin widths of A£ — 140 (Figure |2J and 
A£ = 200 (Figure El see also Table and Table EJ). Com- 
parison of these figures and the actual values and asso- 
ciated uncertainties shows that the cosmological results 
are in excellent agreement. The results are not strongly 
dependent on the bin width, although there is, of course, 
some loss of information at the larger bin width, which 
is reflected in the larger uncertainties. 

We have also tested the effect of varying the residual 
source contribution on parameters. As expected from its 




-0.5 0.5 0.5 10 0.2 0.4 




Fig. 5. — Same as Fig. 12] except for CBIel40+DMR. By compar- 
ing this figure with Fig. |zj we see that the particular choice of bin 
boundaries does not make a significant difference to the parameter 
estimation. This can also be seen from the comparison of the even 
and odd binning cases in Table l4l 




Fig. 6. — Same as Fig.|2] except for the bin size, which is 200 here. 
By comparing this figure with Fig. [5] we see that the bin size does 
not make a large difference to the parameter estimation, excepting 
n b h 2 , for which the Al = 200 binwidth case does not give a useful 
measurement. This can also be seen from the comparison of the 
Al = 140 and Al = 200 cases in Table g] 



small relative contribution, we find very little difference if 
we either assume our standard power value described in 
Papers II and III, allow for a 50% error in that estimate, 
or multiply the standard power by a factor of 2.25. 

We end our discussion of tests of the cosmologi- 
cal parameter extraction from the CBI observations by 
showing projections of the full seven-dimensional like- 
lihood function onto 2D contour plots of the likeli- 
hoods of various combinations of cosmological param- 
eters {weejm, Ofc, Ub, J7a, n s } to illustrate further the con- 
sistency of the various CBI data subsets used throughout 
this work. These 2D contour plots provide a different in- 
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Table 4. CBI Tests and Comparisons 



Priors 




n 3 


n b h 2 




O Ji 2 

J£ cdm"' 


n A 






n b 




h 


Age 


Tc 


CBIol40 
wk-/i 

flat+wk-h 
flat+wk-h+LSS 


i.uu 12 

(1.00) 

(1.00) 


1 - U 'Q.1Q 
i - uo 0.09 


0.0232 
0.0232 
0.024g 


016 

111!! 

00') 


16°° 8 
Q-SK 
U.IOq p. 

u - 11 0.02 


0.43°;| 
0.47°;| 


0.592 
0.542 
0.34° 


22 

M 

li 
12 


0.0832 
0.0682 
0.057g 


053 

\m 
020 


0.58°;" 
0.60°;J| 
0.66«;« 


i3.9|;l 

14.0|;| 


< 0.66 

< 0.65 

< 0.62 



CBIel40 
wk-h 

flat+wk-h 
flat+wk-h+LSS 



°- 96 o 13 L10 o'n 0-016S 
(1.00) I-OSqiq 0.0182 
(1.00) 1.06Hi 0.0202 



8¥8 

007 



0.182 
0.15° 

0.112 



0.372 
0.44° 
0.682 



0.622 
0.562 
0.32g 



I! 0.0592 
f. 0.059° 
g 0.049S 



144 0.60 ' 11 .!•>..,., 

ISi q 5g8:ll ^4 i^'.l 
!?8 0.67§ ; H 1A 



13.3§ } < 0.66 
< 0.66 



14.3^ < 0.63 



CBIol40(£ > 610) 
wk-h 

flat+wk-h 
flat+wk-h+LSS 



1 06°' 15 
i.uo 14 

(1.00) l.lO; 

(1.00) 



1.140- 15 
1 



0.068 



■040 0-047H 
1.05°;" 0.041° 



005 

lis;: 



0.15 

0.18 
0.13 



0.44 
0.41 
0.66 



0.71 
0.62 
0.35 



26 264 - 207 

23 0.188 

H n nso8:03 i 7 ; 

11 U - U °^0.037 



0.59°;" 
0.63P 11 
0.71 



8:11 



13. il l < 0.67 
12.6M < 0.66 
13.l| : | < 0.62 



CBIo200 

wk-h 1.12013 1.14°-12 0.048°-°23 0.220-0° 0.36°- 27 , 0.82° j 2 0.152°-°§| 0.58°-j° 12.8J-0 < 0.67 

flat+wk-/i (1.00) l-07° : ii 0.025° : °if; 0.19° : °f 0.41° : !| 0.61° : |f 0.071° : °| 0.63° : ji 12.6* : | < 0.64 

flat+wk-h+LSS (1.00) 1.04°;J| 0.028°;°jf 0.12°;j| 0.68°;°o 0.33°;jf 0.061°;°i 0.70°;JJ 13.6j;| < 0.59 



CBIdeep 

wk-h 1.09° M 1.16°-|| 0.078°-°|° 0.210-11 . 42 o.3i . 85 o.| .26l84g 0.6ljHg 12.2j-| < 0.67 

flat+wk-h (1.00) 1 - 05 oil °- 050 o'Q34 °- 20 fil4 °- 37 Q ?5 °- 65 fiS - 189 o'll °- 64 Q U 12 - 2 19 < °- 66 

flat+wk-h+LSS (1.00) 1.03[J;lf 0.055§;g|| 0.13[J;g| 0.55g;|| 0.50g;|| 0.187^;^ 0.66g;J| 12.9|;g < 0.65 



DASI+CBIol40 
wk-h 

flat+wk-h 
flat+wk-h+LSS 



1.05°;°| 1.012-11 0.023° 



(1.00) 
(1.00) 



0.992 
l.OOf 



'lit 



0.021 
0.022 



on.] 

881 



0.12° 

0.14° 
0.12° 



0.55° 
0.562 
0.66g 



0.51° 
0.462 
0.33g 



10 0.0778-024 
SI 

?4 

10 "■"''"O.OOQ 



0.0572 
0.0482 



0.56°;i 
0-68«; 



15.21 
13.9° 
13.8° 



< 0.63 

< 0.39 

< 0.42 



DASI+Boom+CBIol40 
wk-/i 

flat+wk-h 
flat+wk-h+LSS 



1.03 



0O.05 

a 0.05 

(1.00) 
(1.00) 



Qt;0-09 

QQ5 
0.940-OB 

0.94 



0.022; 
0.02i; 
0.022: 



0.13; 
0.14| 
0.13 



0.03 

83 
0.02 



0.52 
0.55 
0.63 



0.52 
0.48 



0.37 



8:2? 
8:?4 



0.074°; 
0.058«; 
0.051°°°! 



5fi0.io 
06x8:19 

0.65§$ 



15. OH < 0.52 
14.0° : ^ < 0.31 



1.3.9; 



0.5 



< 0.36 



all-data 

wk-h 1.04°;°1 0.98°;i° 0.023°;°| 0.12°;°3 .55°;17, 0.4»g;« 0.075°;°| 0.56°;" 15.1M < 0.57 

flat+wk-h (1.00) °- 96 Q05 °- 022 Q'oO? °- 13 « 03 - 60 Q'26 °- 42 Q 20 °- 055 Q015 °- 64 0'll 13 - 9 Q 5 < °- 35 

flat+wk-h+LSS (1.00) 0.97[j;g| 0.022g;g^ 0.12g;g| 0.66g;f| 0.35g;lg 0.05lg;gg| 0.66g;g| 13.Sg;| < 0.39 



Note. — Cosmological parameter estimates as in Tabled] except for a variety of data combinations which test and compare results. Only the 
wk-/i, flat+wk-fa and flat+wk-h+LSS priors arc shown. 



sight into the tests of data consistency than do the ID 
plots of the previous subsections, which can be illumi- 
nating. For example, they can show directly whether 
there are isolated multiple peaks in the likelihood surface, 
whereas this information is lost in the marginalization of 
the ID plots. 

In Paper III we showed the power spectra for the indi- 
vidual mosaic fields, and we found that the differences 
are not statistically significant. Here we derive con- 
straints on the cosmological parameters from three sets 
of pairwise-combined mosaic fields. These give stronger 
constraints on parameters than the single fields and thus 
provide a stronger consistency check. We find there is 
good consistency between the three pairs, showing that 
no field is seriously discrepant with the other two, in 



agreement with our finding in Paper III that there are 
no significant differences between the spectra from the 
three fields. The agreement between the three pairs of 
fields can be seen in Figure [7| the results of the pairwise 
splitting of the three mosaic data sets. Displayed here 
are the 1- and 2-a contours for the odd binning of the 
A£ = 140 data. 

The effect of introducing the LSS prior is also shown in 
Figure It primarily reduces the uncertainties in w c dm- 
Note that it is the combination of the CBI observations 
and the LSS prior that accomplishes this. This can be 
seen by comparing Figure [21 in which the uncertainties 
under the LSS prior are large, to Figure [21 in which the 
CBI data have been included and the uncertainties are 
greatly reduced. 
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Fig. 7. — Comparison of likelihood contours for the three pairwise combinations of the three mosaic fields for the CBIol40+ DMRdata. 
Full lines are la and dashed lines are 2a contours. The hatched regions indicate portions excluded by the range of parameters considered 
(see Table EJ. The weak-ft prior applies in the oj c( j m -Cfe panel and fiat+weak-?t prior applies in the other panels. The plots on the left 
do not have the LSS prior, while those on the right do. The pairs of mosaics are as follows: 02h+20h (magenta); 02h+14h (dark blue); 
14h+20h (green). The results for the three fields combined are shown by the dark grey regions (lcr) and the light grey regions (2a). For 
comparison, the results for the alternate binning (CBIel40+ DMR) are shown by the gold contours. Comparison of the plots on the left 
with those on the right show the dramatic impact of the LSS prior on the oj c( j m determination, its role in f^A determination, and its 
relatively small impact in the other variable directions. The dramatic reduction in the uncertainties in u> rc i m is due to the combined effects 
of the CBI data and the LSS prior (see text). (For this comparison, the slightly smaller database in (Lanse et al. 2001) was used rather 
than that in Table [I]) 



A 2D comparison of the parameters extracted under 
the two binning schemes is shown in Figured with the 
effect of adding the LSS prior also shown. These should 
be compared with Figure [2] and Figure [SJ and with the 
values given in Tabled We see here, once again, that 
the consistency between these two binning schemes is 
excellent. 

One of the caveats that is important to bear in mind 
when dealing with 7D-plus parameter spaces is that the 
limits derived from projection onto one or two direc- 
tions in parameter space by full Bayesian marginaliza- 
tion (integration) over the other variables can be mis- 
leading in certain cases. For example, highly likely 
models which exist far from broad likelihood peaks 
may be ruled out. With all of the CMB data, many 
of the variables are well localized and this is not a 



big problem. However, near degeneracies among cos- 
mologic al parameters do ex i st in these inflation model 
spaces fEfst athiou fc Bondl 119991) . The correlations 
can b e quan tified b y considering parameter eigen modes 
(e.g.. iBonl H9M Bond Efstathiou fc Tegmarkl Il997t 
lEfstathiou fc Bondlll999t ILange et al.1120011) which yield 
linear combinations of the C^-database variables that give 
orthogonal errors in the neighborhood of the maximum 
likelihood values. 

For CBIol40+DMR with weak-ft, prior, two combina- 
tions are determined to better than 10% and two oth- 
ers to 15%. The three best-determined values involve a 
predominantly n s combination, with tq and amplitude 
mixed in, and then, somewhat remarkably, a predomi- 
nantly Clk combination, followed by a mix of many pa- 
rameters, including amplitude, Q\ and uj c dm- With all- 
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data, weak-h gives 4 combinations to better than 10%. 
The first two are different mixes of n s and f2fe, next is 
mainly u>b and then a w c dm _ Wb-f2A combination. LSS 
sharpens up the fourth eigenmode especially. 

In conclusion we would emphasize the excellent over- 
all consistency of the cosmological parameters extracted 
from the CBI mosaic and deep observations. It was 
shown in Papers II and III that various subsets of the 
data are consistent with each other to within the levels 
expected, given the uncertainties. Here we have shown 
that the cosmological parameters derived from different 
subsets of the data, and from different binning schemes 
and bin widths, are also self-consistent. We conclude 
that both the CBI data themselves and the cosmologi- 
cal parameter extraction from the CBI observations are 
robust. 

5. OPTIMAL POWER SPECTRA 

In this section, we apply the parameter-estimation 
methods of Appendix to combine power spectra de- 
rived with different bandpowers and with different win- 
dow functions onto a common set of bands. A byproduct 
is the determination of how consistent the power spectra 
are and the values of various experimental parameters 
that can be adjusted to increase agreement (e.g., adjust- 
ments to the flux density scale) . For all of the bandpower 

( s) 

applications, a best-fit model was used for C\ rather 
than flat bandpowers. 

5.1. CBI Optimal Power Spectra 

In addition to the extraction of the optimal power spec- 
trum for the whole CBI data set, we are interested in 
extracting optimal power spectra for subsets of the data. 
Although these optimal spectra are not used in the esti- 
mation of the cosmological parameters, they do provide 
an invaluable means of comparing the various data sets, 
and, of course, of comparing the CBI data to other CMB 
data sets. 

We begin the extraction of optimal spectra by examin- 
ing the pairwise combinations of the three mosaic fields. 
Thus we combine pairs of the three mosaic fields (02 h , 
14 h , and 20 h ). These are separated by about 6 h in RA, 
which is sufficient to treat them as independent fields 
and uncorrelated data sets. When combining the fields 
we do not include separate calibration errors since these 
are common to all fields, and we assume no error on the 
estimate of the residual source component. The NVSS 
sources are projected out as described in Papers II and 
III, and that is included in this treatment. 

The upper panel of Figure ISTT1 compares the combined 
spectra for the three pairs of mosaic fields. The center 
panel shows the combination of the three individual field 
spectra, the spectrum obtained from the joint analysis of 
the three fields (Paper III), and the combination of all 
six mosaic and deep fields. We have separated the upper 
and center panels for clarity, but when plotted together 
they show excellent consistency between all the combi- 
nations we have considered. In the spectra shown here, 
the CBIol40 data have been used and combined onto 
the A£=200 bins. We find excellent agreement between 
the extracted optimal spectra, regardless of whether the 
CBIel40 or CBIo200 data are used. This shows not 
only that the extracted data are self-consistent, but also 



that the band-to-band correlations are being treated cor- 
rectly. We see from the upper panel of Figure 15.11 that 
the agreement of the pairs is robust. This result is re- 
lated to that obtained by comparing the power spectra 
for each individual mosaic field shown in Paper III. These 
show excursions among the power spectra for individual 
fields, although, as shown in Paper III, these excursions 
are not statistically significant. The pair-combined spec- 
tra shown here, with their increased statistical weights, 
indicate again that the excursions are compatible with 
expectations. 

The three deep fields are single pointings with long in- 
tegration times. Since the deep observations have better 
signal-to-noise ratios in the £ > 2000 range, it is use- 
ful to combine the deep spectra with the mosaic spectra, 
which have less cosmic variance at low £. One point to 
note is that two of the deep fields are embedded in the 
mosaic fields. However, the data used for the mosaics 
are only a subset of the data used in the the deep fields 
analysis, corresponding to the typical integration time 
on an individual mosaic pointing, which is about four 
hours. We therefore expect the correlations between the 
deep and mosaic data to be small and we ignore them 
when forming the optimal spectra. The center panel of 
Figure 15.11 compares the combined CBIdeep and joint 
CBIol40 data with those from just the joint CBIol40. 
Both cases use top hat window functions and are mapped 
onto the A£ = 200 mosaic binning. The excess power 
anomaly seen in the deep data is not that evident with 
these relatively narrow bands, and the basic result is 
good agreement between the two. 

The lower panel of Figure 15.11 compares the mosaic, 
deep, and mosaic+deep optimal power using the much 
coarser CBIdeep A£ binning. This is similar to Figure 14 
of Paper III, for which the mosaic data were evaluated 
directly in the CBI-pipcline analysis onto the CBIdeep 
bands. Here we have obtained a similar result by using 
only the CBIol40 data: although there is no detection of 
mosaic power at £ ~ 2200, there is agreement in enhanced 
power levels at £ ~ 2700, although the larger ler mosaic 
confidence region is consistent with a "non-detection". 
The joint deep+mosaic bandpower has slightly smaller 
errors than the deep-only case. 

The combined spectrum shows a clear detection of the 
expected damping of the power out to £ ~ 2000. Thus 
the unique experimental setup of the CBI has further 
validated the cosmological paradigm outlined in § H by 
confirming one of the key ingredients. Note the agree- 
ment in the excess power in the £ ~ 2700 band (though 
errors differ). Further study is needed to confirm this 
excess, and, if confirmed, to determine the source of the 
power (see Paper VI for further discussion). 

5.2. The Overall Optimal Spectrum 

In this section we compare the Cb (band-power) spec- 
trum from the CBI with the spectra obtained from some 
other CMB experiments, and we then combine the CBI 
data with the data from the BOOMERANG, DASI, 
MAXIMA, and VSA experiments in order to obtain a 
new optimal Cb spectrum out to £ — 3500. This repre- 
sents a considerable extension of the optimal spectrum 
beyond the previous limit at £ ~ 1000. 

A comparison of the pre-CBI optimal spectrum with 
the optimal spectrum including the CBI deep and mo- 
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Fig. 8. — Optimal power spectra from CBI observations. Upper panel: Consistency between pairs of CBIol40 fields. Here the CBIol40 
data have been re-binned to the CBIo200 bands for the purpose of comparison with other data sets (see text). At low-4? the errors are 
dominated by sample variance, which is a fraction of the fitted band-power — thus the apparent discrepancy between the pairs in the 
lowest-^ bin, due to the 02h field, is not significant (see Paper III). The best fit model, under the weak-/i and flat+weak-?t priors, to the 



DMR+CBIol40 data (black curve) has {fitot, Q A , tt b h 2 , Q cdm /i 2 



T C }={1.0, 0.6, 0.0225, 0.12, 0.95, 0.025}. The best fit model, under 



the weak-/i, flat+weak-?t and flat+weak-h+LSS priors, for ALL-data (magenta curve) has parameters {1.0,0.5,0.02,0.14,0.925,0}. Middle 
panel: Combined joint 3-field CBIol40 and CBIdeep power spectra denoted by the green triangles; black circles denote the joint 3-field 
CBIol40 spectrum and red squares denote the 3 CBIol40 fields when combined as with the pair cases. The results are stable independent 
of whether even or odd binning is used. Lower panel: The CBIo200 and CBIdeep data are shown together with optimal spectra for deep, 
mosaic and deep+mosaic on the standard CBIdeep bands. 



saic data is shown in Figure As noted above, each 
experiment has unique binnings, band-to-band correla- 
tions, and calibration and beam errors. This makes a 
straightforward visual comparison difficult. To facilitate 
such comparisons, we therefore construct the optimal 
power spectrum for these experiments binned onto the 
A£ — 200 CBI bins. Sample results are shown in Fig- 
ure El along with some best-fit models. We see here that 
the agreement between the different CMB experiments is 
excellent. In the derivation of the optimal spectra shown 
here calibration uncertainties of 10% for BOOMERANG, 
4% for DASI, 3.5% for VSA, and 4% for MAXIMA were 
included, together with beam uncertainties of 14% for 
BOOMERANG and 5% for MAXIMA. For this purpose 
we need to include the calibration uncertainty for CBI, 
which we take to be a conservative 5%. These are all 
incorporated in the spectra, which leads to significant 
correlations among bands associated with the beam un- 
certainty. That is, one must be careful in using optimal 
spectra directly for parameter estimation since, although 



we can compute the Fisher matrix, we do not know the 
likelihood surface in detail. 

We find the optimal power spectrum for 
BOOMERANG, DASI, and CBI requires a decrease in 
the temperature calibration by a factor 1.015 ± 0.04 for 
CBI and 1.05 + 0.05 for BOOMERANG, and an increase 
by a factor 1.01 ± 0.04 for DASI. The data also favor an 
increase in the BOOMERANG beam by 1.07+ 0.04 for 
the odd binning and 1.05 ± 0.04 for the even binning. 
These values are consistent with the quoted errors in 
all cases. Indeed, assuming that the power spectra are 
derived from a single underlying spectrum, which the 
different experiments are sampling in their respective 
regions of the sky, one could use this technique to 
determine the calibrations and beams. For example, we 
note that the determinations for BOOMERANG have 
errors significantly smaller than the quoted uncertainties 
once they have been compared with the two sets of 
interferometer data, which have no beam uncertainties 
and much smaller calibration uncertainties. 
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Fig. 9. — Optimal power spectra: CBI compared to BOOMERANG and DASI. Upper panel: Joint CBIo200 power spectrum compared 
with optimal power spectra using BOOMERANG+DASI data (red triangles) and CBIol40+BOOMERANG±DASI data (dark blue circles) 
in the CBIo200 binning. The optimal spectra correspond to a maximum-likelihood fit to the power in b ands for the three spectra, 
marginalized over beam and calibration uncertainty. The best fit curves shown are the same as in Fig. 15.11 Middle panel: same as the 
upper panel, but using the CBIe200 binning. Lower panel: The optimal spectrum for "all-data", with a finer A£ = 75 binning up to ~ 800, 
and deep binning at large i. The i > 2000 excess found with the CBIdeep data is denoted by the light blue hatched region (95% confidence 
limit) . 



An important caveat to these optimally-combined CVs 
is that power spectra for the individual experiments 
were assumed to be independent. We have already 
commented on this for CBI deep and mosaic combina- 
tions. In addition DASI's fields overlap about 5% of 
the BOOMERANG area, so there is correlation between 
BOOMERANG and DASI which is not taken into ac- 
count in this treatment. The correlation would have to 
be addressed in order to claim absolute accuracy in the 
adjustment of experimental parameters. 

Figure [§| shows very good agreement in the 300 < I < 
1000 overlap region for BOOMERANG and DASI com- 
bined onto CBI points, and with all 3 experiments com- 
bined. (Note that the window functions for the combined 
spectra are top hats and not the Wb£ of the data points 
shown.) We also find the combined spectrum looks quite 
similar when all-data are used. 

The method we used for Figure|5]has also been applied 
to construct optimal bandpowers for all of the above data 
(including CBIol40 and CBIdeep) with finer binning for 
£ < 1000. This variable binning makes use of the high 
quality BOOMERANG data with its intrinsic A£ = 50 
spacing at low £ and it also makes use of the high-£ cov- 



erage of the CBI (out to £ = 3500). The spectrum is 
compared in Paper VI with power spectra computed for 
the Sunyaev-Zeldovich effect, and will not be discussed 
further here. 

6. COSMOLOGICAL PARAMETER ESTIMATES FROM COMBINED 
CMB DATA 

We test the consistency of our parameter estimations 
by comparing with different combinations of CBI and 
other data sets. We also derive estimates for the param- 
eters from the full compilation of data available. 

The consistency in cosmological parameter space at 
the 2a level among these 21 experiments, and the 5 
higher precision ones, is illustrated in Figure 1101 This 
shows that using CBIol40 data with DMR is quite 
comparable at the 2a level with what was achieved by 
BOOMERANG (with the smallest errors), DASI, MAX- 
IMA, and VSA. The combined data, with the LSS prior 
applied as well, gives the bull's-eye determination of Fig- 
ure EH 

Table0]also shows the sequence obtained when we add 
data from CBI's sister interferometry experiment, DASI, 
to CBIol40+DMR. We then add the BOOMERANG 



14 



SIEVERS ET AL. 




Fig. 10. — Comparison of different experiments. 2-cr likelihood contours for the weak-h prior ((^cdm-H^ panel) and flat+weak-h prior for 
the rest, for the following CMB experiments in combination with DMR: CBIel40 (black), BOOMERANG (magenta), DASI (dark blue), 
MAXIMA (green), VSA (red) and "prior-CMB" = BOOMERANG-NA+TOCO+Apr99 data (light blue). Light brown region shows the 
2-<r contour when all of the data are taken together, dark brown shows the 1-tr contour. The LSS prior has not been used in deriving the 
plots on the left, but it has for those on the right. The hatched regions indicate portions excluded by the range of parameters considered 
(see Table HI. This figure shows great consistency as well as providing a current snapshot of the collective CMB data results. Even without 
the LSS prior (or the HST-h or SNla priors), localization of the dark matter density is already occurring, but Q\ still has multiple solutions. 
The inclusion of the SNla and/or the HST-/i priors does not concentrate the bull's-eye determinations much more for the all-data shaded 
case. Note that the expectation of minimal inflation models is that ss 0, n s ss 1 (usually a little less). The Big Bang Nucleosynthesis 
result, wf, = 0.019 ± 0.002 also rests comfortably within the bull's-eye. As expected from the results in Table 151 relaxing the fiat criterion 
has little impact. 



data in order to check for any effects arising from the 
combination of data from completely different experi- 
mental setups. Finally we show the all-data combina- 
tion. 

A full suite of priors for this all-data case is given in 
Table and the corresponding ID likelihood plots are 
shown in Fierure HTI Figure ITSI shows the small difference 
when we use BOOMERANG+DASI+CBIol40+DMR. 
These combined data yield parameters consistent with 
those derived individually: the curvature is close to flat, 
the spectral index is close to unity, and the baryon den- 
sity is near that fa voured by Big Bang nucleosynt hesis, 
L0 h = 0.019 ± 0.002 lIBurles et al.ll200C* lOlive et al.lll999t 
iTvtler et al]l200C|) . 

The tables also show estimates for what we term "de- 



rived" parameters. These are parameters that can be ex- 
pressed as functions of our C^-database parameters. The 
combinations are Q m , Jib, h, and the cosmological age to- 
We calculate the means and variances of these functions 
over the full probability distribution. We have also ap- 
plied the same method to the computation of the statis- 
tics of t-pkj , C p k,j , and ^dipj , Clip, j , and the position and 
amplitude of the jth peak and jth dip in Ci. These are 
discussed in detail in the following section (§ and the 
results from the application to the mosaic data are shown 
in Paper III. We also determine the values of VL m h 1 T,and 
T e ff in a similar fashion, as reported in Paper VI. Paper 
VI also reports on calculations of alternate amplitude pa- 
rameters to Cio, in particular for the as amplitude used 
in the LSS prior. 
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Note. — Cosmological parameter estimates as in Tabic but now for all-data. 
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Fig. 11. — ID projected likelihood functions as in Fig. [2] except 
calculated using all-data (see Table HI . 



The results set out in the tables also show how applying 
the prior restrictions LSS, SNIa, or HST-ft. to the CMB 
data gives compatible results. 




Fig. 12. — ID projected likelihood functions as in Fig.l2l except 
calculated using CBIol40+DMR+DASI+BOOMERANG+VSA 
(the data for BOOMERANG are taken from lRuhl et alJI2002fl . 



7. REHEATING, RECOMBINATION, AND DAMPING 

In this section, we summarize the basic physical effects 
expected to have an impact on the Cg spectra we observe, 
working back from now, through the reionization of the 
universe, § 17.11 and to the critically important regime for 
CBI observations, recombination. We discuss briefly the 
history of computations of the anisotropy in § 17.21 but 
our main goal is to use the parameters of our best-fit cos- 
mological models to evaluate the physical and multipolc 
scales characterizing the decoupling epoch in § 17.31 

7.1. Reheating 
We have used the depth to Thomson scattering from 
a time t\ to the present to, tq = J t ° cdtn e aT, as one 
of our major cosmological parameters, where n e (t) is the 
average electron density, <tt is the Thomson cross section 
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and c is the speed of light. If we assume the Universe has 
been fully ionized below a redshift z rc h, we have 

t c ~ 0.1(c; b /0.02)K l /0.15)- 1 / 2 ((l + z roh )/15) 3/2 . (2) 

A minimum value for z rc h is ~ 5, so tc > 0.03 is ex- 
pected. The visibility to Thomson scattering is defined 
by e~ TC . As long as tq is not too large, Ci is suppressed 
by a factor exp[— 1tq\ on scales smaller than the horizon 
at z le h, and in particular over the regime probed by CBI, 
BOOMERANG, DASI, MAXIMA, VSA, etc. 

The mechanism for reionization is thought to be the 
overlap of H II regions generated by massive stars housed 
in the very small, earliest galaxies to form. Other pos- 
sibilities, e.g., involving particle decays tuned so that tc 
would not be too large, are more exotic and require extra 
parameters for the theory. We know that the stars must 
form, but little about the efficiency of forming the first 
stars. However z re h is necessarily tied to the formation of 
nonlinear gas structures, hence to the power spectrum of 
density fluctuations. For the ACDM models preferred by 
the CMB data, with n s ~ 1, this implies z rc h should not 
be much more than 15, hence our expectation is that tc 
should not be much greater than 0.1. Although the tc 
likelihoods of Figs. llllll2l now fall off nicely beyond 0.2 or 
so, a limit as strong as this still eludes us. The inability 
to determine tc with higher precisi on is attributable in 
part to parameter near-degeneracies (Efstathi ou fc Bondl 
1999). However the fact that we have detected power at 
large £ shows that z ro h cannot be too big: e.g., although 
some pregalactic energy injection at z ~ 50 is still pos- 
sible, it now seems unlikely that it could have led to full 
reionization at such high redshifts. 

The late time reheating described above is almost en- 
tirely a damping effect associated with photons coming 
towards us carrying anisotropy information being scat- 
tered away from our line of sight. There is a small ef- 
fect associated with new scattering sources regenerating 
CMB anisotropies: the reciprocal effect of photons be- 
ing scattering into our line of sight does not appreciably 
add to the anisotropy we observe, so the net effect is 
the damping decline oc e~ TC . The differential visibility 
Vc = de~ TC /din a — e" TC n e aT / 'H (z) , where H(z) is 
the Hubble parameter at redshift z and a = (1 + z)^ 1 
is the expansion factor, defines the dominant regions in 
redshift where the scattering leads to "visible" conse- 
quences. Reionization results in a bump in Vc around 
z rc h with a tail to lower redshift. 

7.2. Anisotropies from the Recombination Epoch 

We now turn to effects associated with the recombi- 
nation of the primordial plasma. The essential ingre- 
dients were w orked out immediately after the discovery 
of th e CMB l|Peebleslll96i IZeldovich. Kurt fc Sunvaevl 
1969). The novel features are the dominant roles played 
by the two-photon decay of the 2s state to the Is state 
and leakage from the Lyman alpha line associated with 
the expansion of the Univer se. Improveme nts including 
helium recombination (e.g., IHu et al.llT995t) and a more 
sophisticated treatment of hydrogen recombination were 
essential for the high precision era we are entering now 
(|Seager. Sasselov. fc Scotll999fl . 

The great simplification afforded by the smallness of 
the primary anisotropies is that linear perturbation the- 
ory can be used and the photon transport equations can 



be decomposed into independent modes characterized by 
a comoving wavenumber fc. Each mode contributes to 
Ci. However, as mentioned in § Q] the transport of the 
photons through recombination involves all of the com- 
plications of radiative transfer as one passes from an op- 
tically thick to an optically thin "Thomson scattering 
atmosphere," compounded by a changing gravitational 
potential. 

Many attempts have been made, in the long his- 
tory of CMB, to deconstruct CMB anisotropies into 
components associated with baryon-photon acoustic 
compressions and rarefactions, the Doppler effect, 
damping, finite decoupling surface width, polarization 
development, and post-decoupling free-streaming prop- 
agation. These have included analytic, semi-analytic, 
and various small- angle and large- angle approaches. 
There were usually two goals: first and foremost to 
understand the physics defining the basic features of 
the spectra; secondly to make quantitative numerical 
estimations appropriate for the computer power of 
the times. Some used ph oton-baryon one-fluid or 
two-flu id app roximati ons fe.g..lDoroshkevich et al.lll988t 
Bondl 119881 119961 iBonometto. Lucchin fc Valdarninil 



Bond 1988, 1996; Bonometto. Lucchin & Valda 
1984: iDoroshkevichl 119881: IHu fc Sudvamal 
Hu fc White! 119961: IPress fc Vishniad I198CI % 



tion 
Bardeen 



Sclia" 

iSilkl IL%1 iStarobinskvl Il988t iWeinber 
others used these methods in conjunc 
with other transport app roximations, (e.g., 



1968 ( unpublished). iBond fc Efstathioil 
1987 : Iboroshkevich. Zeldovich fc bunvacv 



1996t IKaiserl 1985 : 
2002: iPeebles fc Yul p7C : 
12001 albt IWilsoul 11985 : 



1984, 

1978 1 lJungman et al 
Kaplinghat, Knox & Skordis 
VittormfcSil1dlT98¥~ 
Wilson fc SilkHl98l7 

Given the mutual interdependence of the effects, the 
semi-analytic methods can only be taken so far, and nu- 
merical computation of spectra is the preferred method 
for this high precision era of CMB observations. There 
were many groups who developed codes to solve the per- 
turbed Boltzmann-Einstein equations when dark mat- 
ter was present prior to and following shortly after 
the COBE discovery JBondl 119961 IBond fc EfstathioJ 
1984 L ll987tlCrittenden et alJH998lbl: iDodelson fc .Tubas! 
" 19941 lEfstathiou fc Bondl Il98fit iFukugita et alJ | 199(1 
Gorsk L Stompor fc Juszkiewiczll993UGouda et alJl!991t 
Hu et alJll995HKnoxlll995l: iVittorio fc Silklll984 11992ft . 



Most of these solved hierarchies of coupled multipolc 
equations. A speedy, publicly available and widely used 
code for evaluation of anisotropie s in a variety of cos- 
mological scenarios, "CMBfast" (ISeliak fc Zaldarriagal 
1996), using line-of-sight past history integrations of 
CMB anisotropy source terms, has come to dominate the 
scene, used even by those who developed their own codes 
and was checked in detail with a number of these other 
transport codes. It has also had extensions to more cos- 
mological models added by a variety of researchers. An- 
other fast code, "CAMB" ({Lewis. Challinor. fc Lasenbvl 
2000), is based on this technique. A suitably modified 
CMBfa st was used in the c onstruction of the C^-database 
used in lLange et alJ <|2001|) and in this paper. 

7.3. Numerology of Recombination and Damping 

We have seen that analytic, semi-analytic and decon- 
struction attempts have been historically important in 
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quantitative CMB work and continue to be qualitatively 
uselul in understanding how the various effects manifest 
themselves in Ci. Here we shall concentrate on some 
of the scales of relevance for this quali t ative description 
llBcmal Il99fit lEfstathiou fc Bond! 11993 IHu fc SugivanTal 
l1995HHu fc WhitelM 996ft . 

As we take t\ through photon decoupling and recom- 
bination, Tc(ti) grows to extreme opaqueness above the 
rcdshift z dcc where photons decouple, and where recom- 
bination predominantly occurs. The differential visibility 
Vc is sharply peaked for normal recombination and only 
weakly dependent on cosmological parameters. Zd cc is 
defined to be where Vc has a peak. A value ~ 1100 
is obtained for a wide range of cosmologies. For the 
fitot = 1 ACDM cosmologies favored by our CMB data, 
z dcc ~ 1050, with age t dcc « 400000 yr, (comoving) 

— 1/2 

"horizon" scale CTd cc ~ 110 Wm Mpc, and (comoving) 

distance from us Xdoc ~ 0.88 (6000 ui m Mpc). (For Ein- 
stein de Sitter universes, the 0.88 becomes 0.97.) 

The ratio of the density of ordinary matter to relativis- 
tic matter at z dcc would then be ~ 3.3, and the baryon 
to photon density ratio would be ~ 0.8. The sound speed 
at decoupling, c s = (c/ v / 3)[1 + 3 / o_b/(4p 7 )]~ 1 / 2 is lowered 
over the cj \[Z for a pure photon gas because of the inertia 
in the baryons; for w& = 0.02, it is w 0.8c/V3, leading to 

— 1/2 

a sound crossing distance c s Td ec ~ 50u; TO Mpc. Since 
c s varies with time, an appropriate average c s should be 
used, resulting in an adjustment upward of 12%. The 
phase of the waves as they hit the narrow recombination 
band, kc s T dcc , determines the oscillations in Cg , associat- 
ing peak m with a length scale ~ (c s Td ec ) / '{inn). 

To convert the comoving distances at Zd oc to angu- 
lar scales, we divide by Xdcc- The component of a 
wavenumber perpendicular to the decoupling surface, k± 
is associated with a multipole £, where fcj = Xdec^ > 

fa 5.3(1000/^)w„ 1/2 Mpc for our A CDM case. The 
mass in matter enclosed within a perturbation of ra- 
dius r is M = 2.76 x 10 u w m (47r/3)r 3 M . The appro- 
priate top hat radius of a collapsed object that forms 
from waves associated with a band about k is r ~ 2A; -1 
ijBond fc Mversiri996D . This gives a rich cluster mass, 

1.3 x 10 15 w™ 1/2 (1000/£) 3 M Q , for^~ 1500. 

Converting peaks in fc-space into peaks in ^-space is 
obscured by projection effects and the influence of other 
sources such as the Doppler term. The conversion of 
the oscillations into peak locations in Ce gives £ p k, m ~ 
/ m TO7TXdec/ (c s T"dec); where the numerically estimated f m 
factor is « 0.75 for the first peak, approaching unity for 
higher ones. These numbers accord reasonably well with 
the values we obtain when we average over the probabil- 
ity distribution for our C^-database. Using all-data, and 
the flat+wk-Zi+LSS prior, the first five peak locations 
are at 221 ± 2, 536 ± 5, 814 ± 17, 1129 ± 18, 1427 ± 20. 
The interleaving dips are at 413 ± 4, 673 ± 8, 1015 ± 10, 
1310 ± 14, similar to the predicted (m + 1/2) spacing. 
For all-data and the wk-h prior, the positions are very 
similar, but as expected the errors are slightly larger. We 
saw in Paper III that the "model-independent" estima- 
tion of peak positions from the data accord reasonably 
well with these C^-database determinations. 

The electron density falls dramatically through decou- 
pling: the local power law index, p = —d\nn e /d\iia, 



rises from its low and high z asymptotic values of 3 to 
a maximum of about 15, with 10 to 14 typical at Zdcc, 
the range depending upon the model; e.g., for ACDM, 
p w 12. A "Gaussian" width of decoupling in lna/adec 
can be estimated analytically from p: a a ,dcc ~ (p — 1) 
(apart from a small correction factor associated with the 
change of p). For ACDM, a a . dcc w 0.08. (Although the 
distribution is somewhat skewed, a Gaussian fit to Vc 
turns out to be a reasonably good approximation over the 
dominant range, and estimating <j a d ec from the FWHM 
of exactly Vc yields 0.06 to 0.1 in good accord with the 
analytic estimates.) The corresponding comoving length 

scale, Rc,dec ~ a a , dcc [H{z dcc )a dcc }^ 1 , is ~ 7w m 1/2 Mpc 
for ACDM, to be compared with c s r dcc . Because this is 
parallel to the line of sight, i?c,dcc does not project onto 
an angular scale we can observe, but perpendicular com- 
ponents of this characteristic size would have £ ~ 820. 

A combination of viscous (Silk) damping and fuzzincss 
damping diminishes the amplitude of the acoustic peaks. 
The two effects occur simultaneously, intertwined by the 
complexities of the transport, but are estimated differ- 
ently. 

Earlier than decoupling, the photons and baryons are 
so tightly coupled by Thomson scattering that they can 
be treated as a single fluid with sound speed c s , shear 
viscosity (4/(15/,,))/3 7 /(n e o":r), zero bulk viscosity, and 
thermal conductivity k 7 = (4p 7 /(3T 7 ))/(n e crT)- Here 
is 3/4 if Thomson scattering is fully treated and 1 if po- 
larization and the angular dependence of the Thomson 
cross section are ignored. Silk damping has usually been 
estimated using a WKB approximation to these one-fluid 
equations, which results in an overall damping multiplier 
of form e~( k ° DTdBC > I 2 for z > z dcc multiplying the fcc s rd cc 
terms which give the acoustic oscillations. The parame- 
ter <td is an integral of the damping rate involving the 
shear viscosity and thermal conductivity, 



(p-l)(p-l/2)(15/„)(l + ii) 

& 

+ (p-l)(p- 1/2)12(1 + R)*' 



With f v = 3/4 and 10 < p < 14, we have 0.02 < a D < 
0.03 over a wide range of cosmological parameters, 0.023 
for ACDM, with only a weak sensitivity to wj. The first 
term is from shear viscosity, the second is from thermal 
conductivity. For ACDM with ojf, = 0.02, the ratio is 
5 to 1. Although polarization increases an by 10%, we 
cannot determine the damping scale with such accuracy 
with the CBI data. 

— 1/2 

For this an, we get anCT dcc w 2Au> m Mpc, giving 
a scale in ~ 2160. However, the tight coupling equa- 
tions break down as the radiation passes through de- 
coupling, so it is better to treat an as a phenomeno- 
logical factor and match it to numerical results. An 
estimate of the Silk damping scale (given as well by 

iBondl [1996) was fcgjjk ~ 3.8 u> m Mpc, with angular 
size 2.5' and £snk ~ 1390. A more sophisticated phe- 
nomenology of numerical C^-results adopted damping 
envelope functions, exp[— (£/£n) mD ], multiplying "un- 
damped Cis" and provided fitting formulas for £n and 
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m D l)Hu fc Whitd 11996(1 . For the ACDM parameters 
used here, we get £d — 1345, with a power mo — 1.26, 
in good accord with the £snk estimate. Note that the 
falloff is not as steep as the WKB Gaussian would pre- 
dict. When we average £d over the database, we get 
1352 ± 15 for all-data and the flat+wk-Zi+LSS prior. 

The fuzziness damping acts only on fc|| , the component 
of the waves through the decoupling surface: destruc- 
tive interference from both peaks and troughs occurs 
for waves with fc|| Rc,dec > t, but there is none if the 
photons are only received from cither peaks or troughs, 
but not both, the case if oscillations are along the sur- 
face, or if the wavenumbers are small. (The WKB tight- 
coupling solution does in fact calculate a version of fuzzi- 
ness damping along with other transport effects, but the 
fc_i_-fc|| asymmetry is obscured by the truncation of the 
£- hierarchy at such low I: up to Td cc , higher moments 
are strongly damped, but this is not correct as the pho- 
tons pass through Zdcc-) In the Gaussian approximation 
to Vc, fuzziness damping acts on AT/T through a multi- 
plier e _ ^ fc n i?c ' dec ^ / 2 . Because it acts asymmetrically, it is 
not as dramatic a drop as in the WKB case even though 
Rc,dec is bigger than kg-^ k . A simple estimate of a fuzzi- 
ness damping scale angle-averages (fc||i?c,dec) 2 , reducing 
the effective filter to Rc,dec/V3, giving £c,dec ~ 1420, 
similar to ^silk and £d- 

If we use Id = 1345 and the damping envelope to "cor- 
rect" the heights of the peaks and dips determined along 
with their €-space locations, the peak power bounces be- 
tween 4000 and 2500 fiK 2 (except for the first peak), and 
the dip power between 2600 and 1400, "correcting" for a 
significant fraction of the factor of about five raw varia- 
tions. It is certainly an attractive proposition to directly 
translate the CMB data into accurate determination of 
the physical scales operating at decoupling. However, 
the intertwining of transport effects makes the use of a 
parameterized Ce model space a more robust proposition. 

8. CONCLUSIONS 

The CBI provides a unique view of the CMB spectrum 
extending to much higher I than previous experiments 
which have detected primary anisotropics, and well into 
the multipole region of the spectrum dominated by the 
damping of fluctuations at decoupling due to viscosity in 
the photon-baryon fluid and the finite thickness of the 
last scattering region. The CBI observations indicate 
a flat universe with a scale-invariant primordial fluctu- 
ation spectrum consistent with the inflationary model; 
and in addition they indicate a low matter density, a 
baryon fraction consistent with Big Bang Nucleosynthe- 
sis, a non-zero cosmological constant, and a cosmological 
age consistent with the ages of the oldest stars in globular 
clusters. 

These results hold for the whole CBI data set, and 
in addition they hold for a subset of the data restricted 
to 610 < I < 3500. These findings are therefore inde- 
pendent of the spectrum over the ^-range of the first and 
second acoustic peaks, and thus provide independent con- 
firmation of the major results determined by other CMB 
experiments from observations that span the first two 
or three acoustic peaks. This independent confirmation 
of the major results gives much confidence that the key 
assumptions of minimal inflationary models are correct, 



especially the assumption that the primordial fluctua- 
tion spectrum does not have significant fine structure. If 
there were significant fine structure, it is extremely un- 
likely that spectral studies over different ^-ranges would 
yield the same values of key cosmological parameters. 
The good agreement with the results from lower £ also 
demonstrates that the recombination theory of the sim- 
ple model is substantially correct. 

In more detail: the CBI observations from the first 
year of observing, when combined with DMR, gives the 
following key cosmological results, as discussed in 8 JM 
Under the weak-/i+LSS priors we find fl to t = 1-051q'o|, 
and n s = 1.07^;^, consistent with inflationary models; 
O c d m ft 2 = 0.101qq3, and, in addition, identifying the ex- 
cess energy density with the cosmological constant, wc 
find f^A = 0.67^0 i3- When the more restrictive priors, 
flat+weak-/i+LSS, are used, we find £l c d m ^ 2 = 0.11^ ' 02 , 
consistent with large scale structure studies; flf,h 2 = 
0. 024^0 QQg, consistent with Big Bang Nucleosynthesis; 
n m = 6.34 ± 0.12, and fl b = 0.057 ± 0.020, indicating 
a low matter density universe; h — 0.66^ consistent 
with the recent determinations of the Hubble Constant 
based on the recently revised Cepheid period-luminosity 
law; and to = 14.2ti'g Gyr, consistent with cosmologi- 
cal age estimates based on the oldest stars in globular 
clusters. 

These values for key cosmological parameters are in re- 
markably good agreement with those determined in other 
recent CMB experiments. As pointed out above, this 
is highly significant, since the CBI cosmology is based 
on a higher f-range than has previously been used and 
leverages principally off the damping tail region of the 
spectrum rather than the first acoustic peak. 

Another unique aspect of the 500-3500 Grange that 
CBI has probed is that the angular scales correspond to 
the 3D wavenumbers of structures that collapsed to pro- 
duce clusters of galaxies, ranging from those with masses 
as low as 10 14 M Q , to large superclusters with masses 
~ 1O 16 M0. Hence the CBI observations span the whole 
range of masses from small groups of galaxies to large 
superclusters. Thus detection of CMB power in this re- 
gion provides a direct link between the small AT fluctua- 
tions at the time of photon decoupling and the nonlinear 
density amplitudes on those scales today. This provides 
further strong support for the gravitational instability 
picture of structure formation. 

The CBI, BOOMERANG, DASI, MAXIMA, and VSA 
observations are consistent with one another over the en- 
tire range of overlapping coverage in I. The consistency 
between the five data sets, obtained by different experi- 
ments using different observation strategies on different 
parts of the sky, eliminates many sources of systematic 
error as a potential cause for concern. Figure ITU1 shows 
concordance of the different experiments in our minimal- 
inflation parameter space. 

Given the small signal levels that are being studied 
here, it is remarkable that the agreement between the 
recent CMB experiments (TOCO, BOOMERANG, CBI, 
DASI, MAXIMA, and VSA) should be so good. The 
variety of techniques employed, combined with the ex- 
tension of the spectrum to high-i? provided by the CBI 
data and the high ^-resolution data at low multipoles 
provided by the BOOMERANG, DASI, MAXIMA, and 
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VSA experiments, makes for a compelling case that both 
the observations and the cosmological results are robust. 

We have not treated several other parameters that 
could be of relevance to the inflation-based model. In 
place of ACDM, QCDM is receiving much attention, 
with Q an ultra-low mass scalar field, often called 
quintessence, that dominates at late times. Thus £Iq re- 
places Q\ and an effective Q-dynamics is cast in terms of 
a mean pressure-to-density ratio wq — Pq/pq, an effec- 
tive equation of state (EOS). (The dynamics of Q is more 
complex than this, since Q is expected to be spatially as 
well as temporally varying. There is also no good candi- 
date for a theory of Q.) For A, wq = — 1, but wq < —1/3 
would get our patch of the Universe into acceleration. 
wq is not well determined by CMB data and we require 
super nova information to get a useful constraint on it 
fe.g- lBond eralll2000all2002bl) . The CMB by itself also 
is insensitive to the addition of a light massive neutrino 
(HACDM models) since there is only a small effect on Cg. 
LSS can add discriminatory power but it simply shifts 
the result to slightly lower (but still nonzero) ffy (e.g., 
iBond et alJl2000ri iPoeosvan fe Starobinskvlll995|) . 

The influence of a possible gravity-wave component 
will be explored elsewhere but the main result is that 
it is also expected to have little effect on the main cos- 
mological parameters presented here. It certainly would 
have no impact on the angular scales probed by CBI. 

Although our results show a tantalizing drop in the tq 
likelihood beyond 0.2, as described in § 0there is some 
distance to go to get a detection at the ~ 0.1 values 
simple theoretical predictions give for ACDM. Signifi- 
cant early energy injection in the medium cannot have 
occurred. This constrains a number of possible, if not 
probable, scenarios: it is always possible to generate 
many ionizing stars early by amplifying structure forma- 
tion from rare collapses that occurred at high redshift, 
having a non-Gaussian component at small masses, or 
making the primordial spectrum bluer on small scales. 

The dominant feature of the CBI data is the overall 
decline in the power with increasing £, a strong predic- 
tion of the basic theory of photon decoupling, with a 
damping scale now moderately well determined, as de- 



scribed in §0 The way the simplest inflation-based mod- 
els have survived the dramatic extension to higher I over 
what previous experiments probed is rather amazing. As 
well CBI has provided further evidence that the peaks 
and dips associated with acoustic oscillations continue 
to higher £ and are in roughly the right locations given 
by the emerging ACDM concordance model. Thus, the 
CMB results from the CBI and from the BOOMERANG, 
DASI, MAXIMA, and VSA experiments, both in isola- 
tion and particularly in combination, strongly support 
the chief predictions of the inflation paradigm — the ge- 
ometry of the universe is flat, the initial density per- 
turbations are nearly scale- invariant, and the density of 
mass-energy in the universe is dominated by a form other 
than ordinary matter. Simple models in which structure 
formation is driven by topological defects are difficult 
to reconcile with the CMB observations. These conclu- 
sions are considerably strengthened by the inclusion of 
other cosmological results such as measurements of the 
Hubble constant, the amplitude and shape of the power 
spectrum, and the accelerating expansion rate. These all 
point towards a nonzero cosmological constant. 
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APPENDIX 

TESTING THE OFFSET LOGNORMAL APPROXIMATION TO THE LIKELIHOOD 

Whether we are interested in estimating parameters which are cosmological or those defining an optimal spectrum, 
the full likelihood surface is required, not just its Gaussian approximation which is valid only in the immediate 
neighborhood of the maximum. This can in principle only be done by full calculation, though in practice there are 
two a nalytic approximations that have have been show n to fit the one-point distributions quite well in all the cases 
tried l)Bond. Jaffc fe K nox 20 00; iNetterfield et alJl2002|) . We show these work well for the CBI deep and mosaic cases 
too in Figures TA13I and IA14I The simplest and most often used analytic form is the "offset lognormal" distribution. 
This is obtained by taking a Gaussian in the variable z B = \n(q B + q B t ), where the offset q B t is the effective noise in 
the experiment: 

V(q) ex exp[-i - z) B Tf B ,{z zf], 

BB' 

4t = if + + (ai) 

Here F^L/ is the curvature matrix for the bandpowers q B and T^ B , is its local transformation to the z B variables. 
We sometimes use the ensemble average value of FgL, , which is the Fisher matrix, rather than the curvature matrix. 
Other approaches for evaluating q B t are reviewed bv lBond fc Crittenden! |)2001|l Figures \A 1 31 and I AMI show how the 
offset lognormal approximation provides an accurate description of the likelihood in individual bands for both deep 
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Fig. A13. — This shows how well the offset lognormal approximation (long-dashed dark blue) does in reproducing the likelihood functions 
(solid black) for individual bandpowers q B when the rest of them are fixed at their maximum likelihood values. This is for the first 6 of the 
7 CBIdeep bands. The offset logormal approximation with q^ t determined from our CBI-pipeline reproduces the likelihood function quite 
well to beyond 2-a. Gaussian (dotted black) and log-normal (short-dashed red) distributions are shown for comparison, the former a better 
fit in the noise-dominated bands, the latter in the cosmic variance dominat ed bands. The offset log-normal interpolates nicely between 
the two regimes. The equal-variance approximation (dash-dotted magenta) (Bond, Jaffc & Knox 2000) also fits quite well. This result is 
a cornerstone of parameter estimation, whether it be for optimal spectrum combinations of parameterized Ci shapes or for cosmological 
parameter estimations in a C^-database. 



and mosaic data to beyond the la level. It is also evident that the Gaussian and pure lognormal approximations do 
not fit the likelihoods as well, the former working best in low signal-to-noise bands, the latter in high signal-to-noise 
bands. 

To compare a given theory with spectrum CTi{y a ) with the data using equation (|A1| . the model q B 's need to be 
evaluated with a specific choice for ipse'- 

q B = I[C t ip Bt ]/l[C^<p B t], where X\f t ] = £ h^f^ (A2) 

is the discrete "logarithmic integral" of a function fg. The associated bandpower in (/iK) 2 is 

C B = q B C { ° ] (no sum), C { ° ] = Zpj'W] AM ■ (A3) 

As discussed bv lBond. Jaffe fc Knoxl l)2000[K iKnoxl (|1999f) . and lBond fc Crittenden! 1)200 If) , the window function ipst 
used for this purpose (and therefore with a different notation than the ipbt defined above) is somewhat arbitrary. The 
simplest choice is again that of a top hat Xb(fy- We prefer to use the signal-to- noise windows Wb(£) derived from 
the experiment (Paper IV). We have also considered the truncated form Wb{£)Xbi- We find that the cosmological 
parameters we derive in §0]and §[5] are insensitive to which form we use. (For example, the largest change occurs in 
Ub, by 4% in the mean and 10% in the error. This is also evident from Figure I5"T1 which uses ipse — Wbi for the data 
to create an optimal spectrum with ipse = Xbi] the results are almost identical to the original data.) 

If we decrease the width A£ of bins, we develop more correlation between neighboring bands. The offset lognormal 
approximation has only been shown to give a good fit when considering individual bins. We therefore require weak 
band-to-band coupling. We have carried out extensive tests on the dependence of the parameter determinations on 
binning width and positioning. We showed in § 0] that the inferred cosmologies are insensitive to the choice of bin 
boundaries. 



REFERENCES 

Allen, B., Caldwell, R. R., Dodelson, S., Knox, L., Shellard, E. P. S., Bacon, D., Massev. R., Refregie r, A., & Ellis, R. 2002, MNRAS, 
& Stebbins, A. 1997, Phys. Rev. Lett., 79, 2624 submitted I astro-ph/0203134 1. 

Atrio-Barandela, F., Doroshkevich, A. C, & Klypin, A. A. 1991 Bennett, C, et al. 1996, ApJ, 464, LI 

ApJ, 378, 1 



COSMOLOGICAL PARAMETERS FROM THE COSMIC BACKGROUND IMAGER 



21 




c, c, 

Fig. A14. — Same as Fig. IA131 but for a selected spread of CBIol40 bands showing again that the offset lognormal approximation works 
quite well for the CBI mosaic data. Another necessary ingredient is that the bands must not be so closely spaced that the band-to-band 
correlation as estimated by the inverse Fisher matrix is strong since a weak approximation is used for this band-to-band coupling. Both 
the A£ = 140 and A£ = 200 binnings are adequate choices, and deliver similar parameter determinations. 



Bond, J. R. 1988, in The Early Universe (Proc. NATO Summer 

School) ed. W. G. Unruh (Dordrecht: Reidel) 
Bond, J. R. 1996, in Cosmology and Large Scale Structure, Les 

Houches Session LX, ed. R. Schaeffer et al. (Elsevier), 469 
Bond, J. R., & Crittenden, R. G. 2001, in Proc. NATO ASI, 

Structure Formation in the Un iverse, ed. R. G. C rittenden & 

N. G. Turok (Kluwer), in press jastro-ph /0108204 I 
Bond, J. R. & Efstathiou, G. 1984, ApJ, L45 
Bond, J. R. & Efstathiou, G. 1987, MNRAS, 226, 665 
Bond, J. R., Efstathiou, G., & Tegmark, M. 1997, MNRAS, 291, 

L33 

Bo nd, J. R., fc Jaffe, A. 1998, Phil. Trans. R. Soc. London, 357, 57 

(astro-ph/9809043l 
Bond, J. R., Jaffe, A* H., & Knox, L. 1998, Phys. Rev. D, 57, 2117 
Bond, J. R., Jaffe, A. H., & Knox, L. 2000, ApJ, 533, 19 
Bond, J. R. & Myers, S. T. 1996, ApJS, 103, 1 
Bond, J. R., et al. 2000, in Cosmology and Particle Physics, Proc. 

CAPP 2000, AIP conference Proceedings, Volume 555, ed. J. 

Garcia-Bellida, R. Durrer, & M. Shaposhnikov (Washington: 

American Inst, of Physics), 85 I astro-ph/0011379 I 
Bond, J. R., et al. 2000, in XIX Conference on Neutrino Physics 

and Astrophysics, Proc. Neutrino 2000, ed. D. Hallman, J. Law 

& J. Simpson (Oxford: Elsevier) jastro-ph/001138H 
B ond, J. R., et al. 2003, ApT, submitted (Paper VI, 

|astro-ph/ 0205386 I 
Bond, J. R., et al. 2002, Proc. MRST 2002 (AIP), ed., V. Elias, 

R.J. Epp, R.C. Myers, p.15 jasto-ph/0210007l 
Bonometto, S. A., Lucchin, F. & Valdarnini, R. 1984, A&A, 140, 

27 

Brles, S., Nollett, K., Truran, J., & Turner, M. S. 1999, 

Phys. Rev. Lett., 82, 4176 
De Bernardis, P., et al. 2000, Nature, 404, 995 
De Bernardis, P., et al. 2002, ApJ, 564, 559 

Contaldi, C. R., Hindmarsh, M., & Magueijo, J. 1999, 

Phys. Rev. Lett., 82, 679 
Crittenden, R., Bond, J. R., Davis, R. L., Efstathiou, G. & 

Steinhardt, P. J. 1993, Phys. Rev. Lett., 71, 324 
Crittenden, R., Davis, R. & Steinhardt, P. 1993, ApJ, L13 
Dodelson, S., & Jubas, J. 1994, Phys. Rev. Lett., 70, 2224 
Doroshkevich, A. G. 1988, Pis'ma Zh. Eksp. Teor. Fiz., 14, 296 



Doroshkcvich, A. G., Zeldovich, Ya. B., & Sunyaev, R. A. 1978, 

Sov. Astron. 22, 523 
Efstathiou, G. & Bond, J. R. 1986, MNRAS, 218, 103 
Efstathiou, G. & Bond, J. R. 1987, MNRAS, 227, 33P 
Efstathiou, G., & Bond, J. R. 1999, MNRAS, 304, 75 
Freedman, W. L. 2001, ApJ, 553, 47 

Fukugita, M., Sugiyama, N., & Umemura, M. 1990, ApJ, 358, 28 
Gorski, K. M., Stompor, R., & Juszkiewicz, R. 1993, ApJ, 410, 1 
Gouda, N., Sugiyama, N., & Sasaki, M. 1991, Prog. Theor. Phys., 
85, 1023 

Halverson, N. W., et al. 2002, ApJ, 568, 38 
Hanany, S., et al. 2000, ApJ, 545, L5 

H oekstra. H., Yee. H . K. C, & Gladders, M. D. 2002, submitted 

Jasta>ph/0204295l. 
Hu, W., Scott, D., Sugiyama, N., & White, M. 1995, Phys. Rev. D, 

52, 5498 

Hu, W., & Sugiyama, N. 1995, Phys. Rev. D, 51, 2599 

Hu, W., & White, M. 1996, ApJ, 471, 30 

Jaffe, A. H., et al. 2001, Phys. Rev. Lett., 86, 3475 

Jungman, G., Kamionkowski, M., Kosowsky, A. & Spergel, D.N. 

1996, Phys. Rev. Lett., 76, 1007 
Kaiser, N. 1983, MNRAS, 202, 1169 

Kaplinghat, M., K nox, L., & Skordis, C. 2002, preprint 

(astro-ph/0203413) 
Kovac, J. JVL, Leitch, E. M., Pryke, C, Carlstrom, J. E., Halverson, 

N. W., & Holzapfel, W. L. 2002, Nature, 420, 772 
Knox, L. 1995, Phys. Rev. D, 52, 4307 
Knox, L. 1999, Phys. Rev. D, 60, 103516 
Lange, A. E., et al. 2001, Phys. Rev. D, 63, 042001 
Lee, A. T., et al. 2001, ApJ, 561, LI 

Leitch, E. M., Kovac, J. M., Pryke, C, Carlstrom, J. E., Halverson, 
N. W., Holzapfel, W. L., Dragovan, M., Reddall, B., & Sandberg, 
E. S. 2002, Nature, 420, 763 

Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473 

Mason, B. S.. et al. 2003, ApJ, in press (Paper II, 
astro-ph/0205384 I 

Mauskopf, P., et. al. 2000, ApJ, 536, L59 (BOOMERANG-NA) 

Miller, A. D., et al. 1999, ApJ, 524, LI (TOCO) 

Mould, J. R., et al. 2000 ApJ, 529, 786 

Myers, S. T., et al. 2003, ApJ, in press (Paper IV, 
|astro-ph/ 0205385 I 



22 



SIEVERS ET AL. 



N etterfield. C. B. , et al. 2002, ApJ, in press (preprint 

|astro-ph/ 0104460 I 
Olive. K. A.. Steigman, G., & Walker. T. P. 1999, Phys. Hep., 333, 
389 

Peebles, P. J. E. 1968, ApJ, 153, 1 

Peebles, P. J. E., & Yu, J.T. 1970, ApJ, 162, 815 

Peacock, J. A., et al. 2001, Nature, 410, 169 

P earson, T. J., e t al. 2003, ApJ, in press (Paper III, 

astro-ph/0205388 I 
Perlmutter, fcj., et al. 1999, ApJ, 517, 565 
Pogosyan, D., Starobinsky, A., 1995, ApJ, 447, 465. 
Press, W. H & Vishniac, E. T. 1980, ApJ 236, 323 
Pryke, C, et al. 2002, ApJ, 568, 41 

Refregier, A., Rhode s, J., & Groth, E. J. 2002, ApJ, submitted 

|astro-ph/0203131 ). 
Riess, A., et al. 1998, AJ, 116, 1009 

Ruhl, J.E., et al. 2002, ApJ, submitted I astro-ph/0212289 I 
Scott, P.F., et al. 2002, MNRAS, submitted I astro-ph/ 0205380 I 
Seager, S., Sasselov, D. D. & Scott, D. 1999, ApJ, 523, LI 
Seljak, U. 1994, ApJ, 435, 87 

Seljak, U, & Zaldarriaga, M. 1996, ApJ, 469, 437 



Silk, J. 1968, ApJ, 151, 459 

Starobinsky, AA. 1988, Pis'ma Zh. Eksp. T eor. Fiz. 14, 394 
Szalay, A. S. et al. 2001, submitted (preprint astro-ph/0107419 I 
Turok, N., Pen, U., & Seljak, U. 1998, Phys. Rev. D, 58, 3506 
Tytler, D., O'Meara, J. M., Suzuki, N., & Lubin, D. 2000, Physica 

Scripta, submitted I astro-ph/0001318 I 
Van Waerbeke, L., Mellier, Y ., Fello, k', P en, U.-L., McCrac ken, 

H. J., & Jain, B. 2002, A&A, submitted ( astro-ph/0202503 I. 
Vittorio, N., & Silk, J. 1984, ApJ, 285, L39 
Vittorio, N., & J. Silk, J. 1992, ApJ, 385, L9 
Weinberg, S. 1971, ApJ, 168, 175 

Weinberg, S. 2001, Phys. Rev. D, 63, 123511 ( astro-ph/0103279 I 
Weinberg, S. 2001, Phys. Rev. D, 64, 123512 ( astro-ph/0103281 1 
White, M., Scott, D., & Silk, J. 1994, ARA&A, 32, 319 
Wilson, M. L. 1983, ApJ, 273, 2 
Wilson, M. L., & Silk, J. 1981, ApJ, 243, 14 

Zeldovich, Ya. B., Kurt, V. G. & Sunyaev, R. A. 1969, Sov. Phys. 
JETP, 28, 146 



